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THE ABSOLUTE STABILITIES OF CERTAIN CONTROLLED SYSTEMS 
Yu.S. Sobolev 


(Leningrad) 


A problem [1,2] of comparison in the theory of absolute stability for automatic 
control systems has been solved. Some examples are given. 


1, A problem, which is formulated in a more general form below, has been posed in [1, 2]. 
We have two sets of differential equations 


Let these sets satisfy these conditions: 


a) solutions to (1) and (2) exist and are unique in the region of every point throughout the space of vari- 
ables xy,..., Xp; 


b) the point (0,... , 0) is a unique state of equilibrium for both systems; 
c) the function gy, belongs to class A [1]: 


(2x) > 0, > 0, when and 9, (0)=0 (k=n—m+1,...,n) 
d) (Xy Xn) = 0 only when 


e) by and hy are real numbers, and hy € (0, +0); 


f) the square of the “distance” of (2)'s image point from the origin dies away strictly monotonically for 
any hy €(0, +00), and when t + +o.* 


Prove that the square of the “distance” of (1)'s image point from the origin also dies away strictly mono- 
tonically. 
Or or It 
The proof is as follows. We multiply the {-th (k-th) equation in (1) by r aA? ia) and sum the result- 
ing equations from 1 to n: 


*By “distance” (henceforth denoted by r) we mean a positively defined continuously differentiable function of 
n variables, 


1d 


k=n—m+1 
We do the same with (2): 


Consider (4). Condition "f" implies that the expression on the right is negative for any values for which 
hy, € (0, +@) (k= +1,... , throughout the space of the variables xy, . . . , Xp. 


We use this to find the signs of the terms on the right in (4). We assume that all the hy (k=n—m +1, 
- » M) are sufficiently small, and so we have that 


Dir <O* (5) 


for any Xn. 


Further, we make one of the h, e.g., hs, large, and all the othershp(p=n-m+1,...,s—1,8+1, 
. , n) small; we have 


Now s is arbitrary, and so we can say that 
or 
+» tn) <0 (k= n—m+1,...,n) (7) 


for any Xy,...., Xn- 
Consider now the product 


(a4, +, Pe [on (24, +++ 5 (8) 


Condition “c" shows that this is always positive. 


Henoe the slgn of - + » Xn)] throughout the space of x,, . . . , Xn is the same as that 
or 


of ++,» Xp). From (6) we have that 


*In fact, for any set of fixed values xy*, . . . , Xp* we can always either put down a set of values hy,*(k = n — 
—m+i1,..., Mm) such that 


kean—m+1 
or can readily show that 


or 
= 0, 
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for any Xy,..., Xn. 


Now from (3), (5) and (9) we get 


ar 
=1 k=n—m-+1 


Hence the square of the distance of (1)'s image point from the origin dies away strictly monotonically 
throughout the space of the variables xy, ... , Xp. 


Note. In the space xy, ... , Xm we consider a set of differential equations 


>) 9457; (i=1,...,m), = const, (10) 


(11) 


and a form derived from (10) by multiplying the i-th equation by roe , and summing the resulting expressions; 


1 d(r%) ar 
i j=1 
where r? = > bj;2, 18 a definitely positive quadratic form. The coefficients in this form can be so chosen 


that the cummiens derivative of f,with respect to t, will, from (10), be negative when t > 0 if, and only if, {|All 
has eigenvalues with negative real parts. 


Now we have assumed that r is an arbitrary function; hence we use the above proof and note to justify our 
assuming that there is a function * whose complete derivative is, on account of (1), negatively defined, by virtue 
of the Rous-Gurvits conditions. 


2. Consider some examples. We assume that the local conditions of existence and uniqueness for the solu - 
tion are complied with, and that the point (0, . . . , 0) is a unique state of equilibrium. 


Example 1. Find the absolute stability bound for a control system with one controllers; 


+ Vin =— Nya, 
Tits + = f (0), 
where o = ayy + + — and f(a) is of class A [1]. 


(13) 


The proof above shows that the absolute stability region for (13) is the same as the asymptotic stability re~ 
gion for the system 


+ Un — Nya, 
T + = h (ay: + + — (14) 


where h€ (0, +00). 
We use a method given in [1] to find the bounds to the asymptotic stability region. 


*Such systems have been dealt with in [1, 4-6] and elsewhere. 


n 
j=1 
ain+**?nn 
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The characteristic equation to (14) is of the form 


+ + + (WU + hrT2 + hNG*) + (rUR4-NEh)d +h Na = 0. 
If U, W, N, E, a, r > 0, the following inequality is the sole important one: 


(hrU + hNE)((TPU + (WU + + hNG*) — (hrU + hNE)| — 
— (T?U + WT?)*hNa > 0. 
The equation for the bounds to the asymptotic stability region for (14) is 
+ WTS) + NG*) (rU + NE) —(rU + + 
+ + (rU + NE) WU —(T?2U + WT2)*Na] = 0. 
The envelope to the family of curves described by h is given by 


T2W (NG? + rT?) = NT? — G*U), 
(15) 


(rU + NE) WU = (TU + WT?) Na. 
Hence (15) is the bound to the absolute stability region for (13). 
Example 2. Find the absolute stability bound for a control system with one controller [1]: 


+ + ky: = — 
Wi + Sy2 = f (2), (16) 


where f(o) is a function of class A, and o = ayy + Eyy + @¥y— 12. 
As before, the absolute stability region for (16) is the same as the asymptotic stability region for the system 


+ Uj + = — Nya, 


Win + = h (ay: + + Gy — rya), (17) 


where h € (0, +00). 
The characteristic equation to (17) is of the form 


+ [WU + (S + hr) + + [Wh + U (S + hr) + 
+ NhE]> + [Nah + k(S + hr)] = 0. 


If U, W, S, N, E, k, a, r > 0, the following inequality is the sole important one: 
[WU + 72 (S + hr) + NG*h] [Wk + U (S + hr) + NhE] —[Nah + k(S + hr)| WT? > 0. 
The equation for the bounds to the asymptotic stability region for (17) is 


h? {(T2r 4+ NG*) (rU + NE)| +h{(WU + T2S) (rU + NE) + (kW + US) (T2r + NG) — 
— WT? (Na + + (WU + T2S) (Wk + US) —WT2kS = 0. 


The envelope to the family of curves described by h is given hy 


[(rU + NE)(WU +T72S) + SU + NG*) + NW (G% — = 
= 4 (rT? + NG) (rU + EN)[SU (WU + 72S) + kUW?). (18) 


The last equation is the bound to the absolute stability region for (16).* 
*Other examples of such nonlinear control systems are dealt with in [1] and [5]. 


the 


the Sc 


*Simi 


wh 
syste 
harn 
| 
= 
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Example 3. Consider the system 


= (x) + 
= + (19) 
where f,(x) and f,(y) are functions of class A. 
The absolute stability region for (19) is the same as the asymptotic stability region for the system: 


Y = + aeahay, (20) 
where hy€ (0, +0), hy € (0, +00). 
The characteristic equation to (20) is of the form 


— + A + — 491413) = 0. 


The zero solution to (20) is asymptotically stable if 


Agy < 0, Age < 0, < 0. 
Hence the zero solution to (19) is absolutely stable, because these inequalities are complied with. 


SUMMARY 


The results here show that normal analysis methods for linear systems can be used with many nonlinear 
systems, even if there are several nonlinearities. Popov [9] has observed that these results should be of yalue in 
harmonic linearization applications. 


The author wishes to thank A.M. Letov and B.V. Shirokorad for their assistance in compiling this paper. 
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AN INTEGRAL ESTIMATE FOR SELECTING THE OPTIMUM PARAMETERS 
OF AN AUTOMATIC CONTROL SYSTEM WITH A GIVEN OVERSHOOT 


A.Il. Tupitsyn 
(Khar'kov) 


A combined integral estimate is used to find the optimum parameters of an auto- 
matic control system with a given overshoot. Formulas are given which define the re- 
sponse time approximately. Examples are given. 


If the integral-square estimate J: = \14@ |'dt, where A(t) is deviation of the controlled quantity from 
0 


a new steady value, is used to choose the optimum parameters of an automatic control system with a given over- 
shoot, it is found [1-3] that the transient response is oscillatory and shows a large overshoot. That this is bound 
to be so has been demonstrated [4]. Large overshoots are sometimes impermissible or undesirable, and so other, 
more complex integral estimates are sought to give systems with better damping. A combined estimate (test) 

of this type is (1, 2, 5) 


Jy \ [A (t)]2dt + 
0 0 


. The extremum in this test is an exponential function of time-constant T, so we may expect to get a highly 
damped system if the parameters are chosen in accordance with this test. It is, in practice, difficult to use J,, 
however, because we have to determine T in every separate case. 


Here we deal with a combined test, Jy; = Jy + cc /ASJiJo3, which gives optimum parameters for a given over- 
shoot. Here a is a dimensionless parameter that defines the extremal transient, in particular, the overshoot. 


We can find how the overshoot depends on a by considering the turning points in Jy. 
The argument that gives us Jy is as follows: J, has as extremum the function A(t) = Age” /T, in which T 
can be found from J: 


T= 2J,/M4 


We substitute for T in J, and get 


Jn = 
0 


It will be seen later on that this is, at its turning point, a function which is nearly exponential in form. 
The test can be generalized, however, if we replace the 4 by the dimensionless variable a. 


We then have 


Ju +2 (1) 
AS 


If now we are able to show how the transient response is related to a by considering the extrema in Jy, we 
may say that J, can serve to choose the optimum parameters for a given overshoot. 


The parameters given by Jy; are optimal in the sense that the transient is briefest, because J, is proportional 
to the time the transient lasts, for given overshoot and Ag; hence, minimum Jy, means minimum response time. 


This method, with J,), is not the only one possible. Cvher tests can be derived by the method of dimensions, 
combined with giving dimensionless forms to the parameters that define the shapes of the extrema. 


J, has dimensions [Af [t], and Joy has 

Consider J; + alt Addex: we determine y and z from dimensional principles: 
({A]? (Ay [A]? = [A]? 

Hence, 2y +z+2=2,y—-1=1, soy = 2, 2 = 4, which gives Jy; above. 


The Theoretical Extrema in Jy, = Jy + (a/A$)J3q1 


The theoretical extrema in 


Ju 
0 


are given (see Appendix) by (17) and (20), in which B = f(a) and B, = f(a) are given by (19) and (21). 
Figure 1 shows 8 = f(a). Here a varies from 0 to o as B goes from 0 to 2,176. The minimum in J, oc- 


curs in this range. r 
SHS 
a 
| 
Fig. 1. The curve B = f(a). Fig. 2. The curve By = f(a). 


| 
i 


Now a varies from —o to —4 as 6 varies from 2.176 to a, and we get the absolute minimum in Jy. 
Figure 2 shows 6; = f(a). The relation is double-valued. The minimum in Jy corresponds to By from 0 to 0.76, 
whereas values between 0.76 and 1 give a relative maximum. 

Figure 3 shows 


1 


for several values of 6 and By. 


Trapezoidal frequency characteristics [6] give us A(r)/Ao, which is shown in Fig. 4 (Fy(v) for B = 4.4 ts 
not shown in Fig. 3). 


20} | 
is 
10 
hte) 


Fig. 4. The curve A(r)/Ag for several values of a. 


Figure 5 has been drawn up from Fig. 4, and gives us a as a function of overshoot o. 


The response time can be found roughly when the parameters corresponding to the minimum in Jy; have 
been found. 


. Fig. 3. The curves Fy(v) and F,(v) for several values of 8 and By. 
| 
a=0 


The time r, from Fig. 4 is used to find wo; t, = 


= Tr/We. 
zs NoW We is given by 
te = BQ = (2) 
where Q¢ is found from (13) (see Appendix): 


4JuBarctg*® _ Ar th* By 


4 
- a 


Fig. 5. Overshoot as a function of a. We (0) are tg (0) Av th py’ (5) 


A?nB*A An Bt A, 8 
2711 (1 + B*) arc tg* B (A — 2B) (1 — B2) Ar th B; (A;—2B,) (6) 


Here Tu denotes the minimum and Ty Jog and Fe(0) denote the values of Jy, Joy and F,(0) that corres- 
pond to Jy. 


Here (3)-(5) derive from (18) and (16) (see Appendix). We have derived (6) by substituting for J, Jy, and 
a [from (18) and (19) in the Appendix] in (1). 


All these expressions ought, in theory, to give the same we. The real extrema differ from the theoretical 
ones, however, and so (2)-(6) will give differing values. The mean we should be used, 


The Estimate Jy, = Jilos 


The relative importance of J, in Jyy decreases as a becomes large. We may expect that the transient re- 
sponse corresponding to 8 = 2,176 will be given by Jy» = Fidot- 


In fact, we find that the analogous variation problem for Jy gives us (16) and (17) again (see Appendix), 
with 6 = 2.176, Then Qg is given by 


Qe (1) 
We again use (2)-(4) and (5) with 6 = 2.176 to get we. Now (6) is replaced by 
0.84 Af (8) 
64 arc tg® Jus 


Examples of Uses of Jy, and Jy, | 
Example 1. Consider a second-order system. Let the Laplace transform of the output be 


pP+2kp+i p 


Then Jy, is 


44) 


is 4 
0) 
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We calculate dJ,,/dk and equate it to zero, and hence get 


2+a+V (2+ a)? + 3a(4+4+ a) 


Hence I, 1s real and positive for all positive values of a, and for a between —a and —4 (which agrees 
exactly with theory). The small negative values of a stretch from 0 to —0.27 (the real result differs somewhat 
from the theoretical). In this range , actually has two positive real values, of which the larger corresponds 
to the minimum in Jy, and the smaller to the relative maximum. 


We get J» from 


(4k* +4)? 
Ju= 


The minimum occurs with kop; = 0.867. 
This value of Kopy follows from the general formula for the limit when a - o. 


Figure 6 shows extremal curves for A(t) = 1 — x(t) for different values of a. 


A(t) 
10 


ag 


a7 


jes | 
a 
\ 
\ 


a 


x=05 
a-0 


Fig. 6. A(t) for example 1. 


We now compare the actual and theoretical response times for Kop = 0.867. 
We find we from (2)-(5) and (8) as 1.54, 1.73, 1.37 and 2.16, respectively. The mean we = 1.71. 


Then, from Fig. 4, we find T.. 9%, the time at which A(r) deviates by not more than 10% from Ag: 
T ss 10% = 4.64. Hence, the theoretical time is 


19. 


Thi 


then 


sp 
the 
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a 


005 sec 


Fig. 7. Jy, = f(T) for example 2. Fig. 8. x(t) for example 2. 


tss 10% = 4.64/1.71 = 2.7 sec. 


The actual time given by Fig. 6 is 3.25 sec. 


Example 2. We consider the optimal parameters for the stabilizing transformer in a self-adjusting motor - 
speed control. The circuit, the transfer function, the numerical values of the main parameters, and the expres- 
sion for Jy, are given in [4]. 


For Jo, we have 


We take a = 4 in Jy, in order to get only a small overshoot. 


Figure 7 shows Jy; = f(T), which implies that J,, is a minimum when T, = 0.08 sec. Figure 8 shows the 
corresponding x(t) = Xss — A(t) curve. The control time, t; 19%, is 0.21 sec. 


We now calculate the theoretical time. We have J, = 0.073, Joy = 3.76, Jy, = 0.2. The values of we are: 
19.0; 19.8; 17.8; 18,9; 21.6. The mean is we, = 19.4. 


The theoretical time is t, = 3.88/19.4 = 0.20 sec. 


SUMMARY 


1, The integral estimate Jy; = J, + (a /A4)F4Jox is examined, and it is shown that this estimate can give op- 
timum parameters corresponding to a preset overshoot. 


2. Jy, and Jy are used to find the optimum parameters for linear stationary systems; the actual results for 
the overshoot and response time agree well with those found from the theory. 


The estimates can also be used in statistical dynamics problems. 


APPENDIX 
The Theoretical Extrema of = Jy + (a / AS) 01 
We find the theoretical extrema of J4; for various values of « by the method given in [4}. 
Let A(t) have a cosine Fourier transform [7] 


=A (1) coset dt, (9) 
0 


then A(t) is given by 


| 


2 =] 
F, (@) cos deo. (10) 
0 
We differentiate (10) with respect to t, and get 
(t) = F, (@) wsin teo deo, 
0 
We use Parseyal's relations for the Fourier transform, and get 


We pass to the final limits of integration with respect to w, and use the approximate expressions for the in- 
tegrals given by the formula for the mean rectangles, and get Jy; a8 


The extra condition 
2 


implied by (10), used with Lagrange multipliers, gives us the function 


where A is a Lagrange multiplier. 
We find dF‘/dF,; and equate it to zero to get 


The continuous function F, (w) then is 


F,(@) =— 


where J, and J, denote the values of J; and Jgy that correspond to the extremum in Jy). 


We use the notation 
+ 2a 
(14) 
and the new dimensionless variable 
v=—, (15) 
c 
to get the expression 


whe 


Then (12) gives us A in the form 
T+ 
hence 
Ap are tg B 
Thus F,.(v) is 
F,(v) = (16) 
Then (10) gives us A(r) as 
1 
(17) 


Hence (17) gives the extremum in J; as a function of the dimensionless variable r = wet. 
We relate 8 to « by finding expressions for J; and Jos from (11) and (16) 


280A 
Ge, arc tg? B = | 4B are tg? B 


Amarctgh + B= arctgp — 


Substituting in (13), we have 


A‘ 420 A‘n*AB 
arc 


16, arc tg*B 


whence, from (14), we have 
16 arc tg* B 
(19) 


This variational problem also has solutions if 6 is imaginary. 
We put B = jB, in (17) and (19), and get 


1 
cos tvdv 0 
0) = \ (20) 
16Ar th* B 
(A; 1)’ (21) 


where A, = Arth B, + 8;/(1 — 84), By = Arth 8; — B,/(1 — 89), and B, is a real quantity. 
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USE OF AN INTEGRAL-SQUARE ESTIMATE 
TO DETERMINE THE OPTIMAL PARAMETERS 
OF AN AUTOPILOT WITH RATE FEEDBACK* 


V.D. Matytsin and V.A. Ryapeloyv 
(Moscow) 


The stabilization of an aircraft by an autopilot with rate feedback is studied. 
Analytic formulas are deduced which define the optimum values of the transfer num- 
bers for the autopilot in terms of the aircraft's aerodynamic parameters and of the 
flight conditions. Lag in the controls is allowed for. The integral-square estimate is 
minimized to yield the optimum stabilization parameters. _ 


To define the optimum parameters for the autopilot is one of the most common tasks in designing aircraft 
controls, Aircraft have large speed ranges, and their above-ground heights vary greatly; the calculations needed 
can therefore be lengthy and troublesome. 


An analytic method of defining the parameters (i.e., the transfer numbers) is thus of some interest. 


The problem is as follows; the aircraft's aerodynamic and geometrical parameters are known (e.g., for 
stabilization in a longitudinal plane), and so are the flight conditions; and from these data the optimum transfer 
numbers must be found. The work involved should be minimized. 


This analytic problem is solved here for an autopilot with rate feedback. The idealized law is of the form 
TAS + = k,Ad + k,Ad + k,A5, 


where AS = 9 — 9p is the deviation of the course angle from the set value, Aé is a variation in the rudder setting, 
T specifies the lag in the control, and kg, ky and k, are the transfer numbers. 
The motion about the center of gravity, when the general motion is longitudinal, is described by 
AS Ab+ - Aa + = — hy Ad, (1) 


Here Ac is a variation in the angle of attack, 


4 
2) 


*Read at the seminar on automatic control at the Institute of Automation and Remote Control in 1958. 


y2 
ms 
(4) 
Sb, 
(5) 


In (2)-(5) p is the density of the air, V is the speed, S is the wing area, b, is the mean aerodynamic chord 
of a wing, Jz is the aircraft's moment of inertia, m¥/* is a damping coefficient, me is the aerodynamic damping 
coefficient, m is the static stability coefficient, and mé is the aileron coefficient. 


We assume that the direction of the general motion does not change while the transient about the center of 
gravity is dying away, and so we can assume that 


AS~Aa and 
We then have 


where AS + AD + = — 


__ A fourth-order equation describes the closed-loop system formed by aircraft and autopilot: 


TASY + (1+ +(q— + The) + (he + AY + (6) 


Gurvits' determinants for the system are 


bo = (1+ + take + Tha) (he + hike) —T (ha + — (t+ 7), 
A, + + The) —T (Ita + kyke). (8) 


We take the minimum in the integral of the deviation squared (the integral-square deviation) as criterion 
for the optimal choice of transfer numbers (the integral test). 


This indirect method may give transfer numbers that imply excessive oscillation, but which are suitable 
for subsequent electronic analog studies [1]. 
The integral-square test is Jo =\A0*dt , and any of the usual methods [1-3] can be used to derive it. 
0 


We assume that we have at t = 0 a typical initial perturbation A® = A», AS = AS = 0, and that as t + 
the stabilization process goes to completion; the integral test then takes the form 


+The), (14 (ik, + ) Tay (1 + 
Wag 


(9) 
(ky + kiks) (1+7-) —T (k, + kikg)® (14+ 


J 
+ 


Then (7) and (8) give 


(10) 


Hence, (10) gives us the integral test in terms of the autopilot’s transfer numbers, and of the aircraft's aero- 
dynamic parameters. 


This rather cumbersome form can be simplified by applying the first condition for the test to be a minimum: 


This condition is 
A, [40 — koks (1 |= 0. (11) 
Now, if Gurvits’ A; is zero, Jy has a maximum (the system is not stable; the integral diverges). 
The test will have a minimum if 
= kok, (1 + 12 
= kekts(1+ (12) 
Then (12) gives 
T 
k, + hy (1+ 


We use (12) to simplify the form of the test: 


From (14) we readily get the second condition for the test to be a minimum 0J,/k, = 0: 


ka + kyks = 2 ( 1+ Te) 


or 
% 
Ta) 


k, ky 
Jo + + 
) 
To 
405 


We use (15) to get ky: 


T 
1+ 


Then, with ky and ly given by (15) and (16), the expression for J, becomes 


T 


Thus k, and k, can be found in terms of ky and of the aircraft's aerodynamic parameters. It remains to 
find an analytic condition to guide us to a reasonable choice of kg. 


Now (17) shows that the test follows a cubic hyperbola law (Fig. 1) with respect 
1 toe 


The test varies very slowly with kp in region AB, and so there is no point in in- 
creasing kp further. Now the drive power must be increased in order to increase ky 
(the rate coefficient for the rudder drive), while the transient response is not much af- 
fected. 


0 i The J, = f (kg) curve has a characteristic point of maximum curvature. The best 
Fig. 1. ke is somewhat greater than the k, that corresponds to this point of maximum curvature 
on Jy = f(y). 


The curvature is given by 


The maximum curvature occurs at a ky given by 


i+ 
komax = (—*2) = 0.9 (18) 
The criterion for ky is the inequality 
1+ 
(72) . 


Our analog studies have shown that such a choice of ky is quite satisfactory. 
Hence, the optimum parameters of an astatic autopilot described by the equation 


TAS + = k,Ad + k, AS + 
are given us by expressions 
for Kg 


We 


T 
To 
: 
» 
(2) ; 
by 
406 


= 


we 
= 


qa 


for ky 
ky = 0 
and for ke 
Qa 


The values are refined by using an electronic analog before design values are suggested. 


Consider, for example, the optimum transfer numbers for stabilizing a hypothetical aircraft in banking 
We assume that ky = 0, kg = 30 sec™, 1/Tp = 2.5 sec”! and T = 0.1 sec. 


Figure 2 shows how the integral test varies with ky. Then kgmax is 0.5 sec™* 


We make the transients tolerable by increasing kamax fourfold, and get ky = 2 sec™", We then use (15) to 
get the optimum value for the derivative coefficient, k, = 1.1. 


Figure 3 shows how the normalized integral test varies with k,. 

We use (16) to find k, (the transfer number for the second derivative of the control signal) as 0.16 sec. 
Figure 4 shows how Jy varies with ly. 

Figures 5 and 6 show transient responses for ky = 2, with various combinations of values for ky and ky. 
It is clear that these formulas give transfer numbers that correspond to reasonable control response. 


Some results given in [5] may serye as the second example illustrating our method. An autopilot described 
by 5 = keAS + kiAS + lg AS was considered, and it was suggested that a reasonable method of constructing sets of 
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integral-square deviations numerically should be used to define the transfer numbers. Many calculations are 
given in this paper; these show that the integral test has a sharp minimum when expressed in terms of k, (Fig. 7), 
and that J, varies hyperbolically with ky (Fig. 8). 


The aerodynamic coefficients are not given in the paper. 


SUMMARY 


1, Approximate formulas are given for the transfer numbers of an autopilot with rate feedback; these for- 
mulas are in terms of the flight conditions and aerodynamic parameters: 


1+ 
T | k, 
F .. 
/s 
ky = ko 2rky 
8 


2. The formulas can be used to give good approximations, which may be used in making a final choice for 
the numbers when the complete set of equations is examined on analogs. 


The authors wish to thank A.M. Letov and B.V. Shirokorad for very valuable advice. 
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THE DYNAMICS OF A VIBRATION LOOP 
IN AN ELECTRICAL SERVO DURING FREE OSCILLATIONS 


IL.N. Krutova 


(Moscow) 


The vibration loop has two parallel control channels, which have a common 
feedback. Three first-order equations describe the motion, The motion is studied in 
the phase plane, by point-transform methods. Bifurcation relations are found for the 
parameters; these relations define the bounds to the various forms of motion in the 
parameter space, 


Vibration servos are often used in automatic controls. The vibration loop is an important part of any such 
servo, because it effects vibrational linearization. Internal loops are often used for the purpose; these loops ef- 
fect local feedback between several amplifying units. The feedbacks have lags. 


The rate of change at the servo's output is usually very much less than the rates of change within the vi- 
bration loop. The system can thus be divided,in principle, into two independent parts; the vibration loop is con- 
sidered as a closed dynamic system, and the seryo as a whole is treated as though the vibration loop were a static 
sampled-data unit. 


The stability of any such loop must be examined. Such a loop is essentially nonlinear, and so the stability 
problem is one in nonlinear theory. Andronov's [1,2] point-transforms method, and phase-space methods, are 
used here to solve the problem. Several papers [3,4] deal with phase-space methods used to solve nonlinear vi- 
bration-loop problems, but in all cases the loops had unbranched input circuits. 


The Equations of Motion 


We consider the vibration loop in an electrical autopilot [5]. Figure 1 shows the loop's electrical system. 
The input U is fed to the grid of tube 1, which amplifies it and feeds the dovble triode 2. The plates of triode 2 
are ac-fed in opposite phases, The tubes conduct during the positive half-cycles only. The second half of the 
tube 2 conducts when the signal changes in sign. Relays P, and P,, with capacitors Cy and C,, smooth the current 
pulses in the plate circuits. The relay contacts kP, and kP, operate the effector, and close the feedback circuit, 
whose signal is fed to the grid of tube 1, where it adds to the input. The circuit Ry Cy in the feedback loop 
causes a lag in the feedback signal when the relays operate and release. 


The following assumptions are made in solving the mathematical problem for the system. 
1, That the first stage is linear and lag-free. 


2. That the second stage is lag-free, with static characteristics F*,(Uy) and F%,(Uy) for the two halves of 
tube 2. Figure 2 shows F*,(Uy) and F*, (Uj) (U, is the voltage applied to the grid of tube 2). 


3. The plate circuits of tube 2 have inductance and resistance only, with different values for positive and 
negative signals, i.e., when tube 2 is conducting and not conducting. 


4. ‘The contacts’ effects are described by the nonlinear functions Up (ig) and Uy (ig) shown in Fig. 3 (1, 
and i, are the currents in the plate circuits of tube 2). 


4 
To sensitivity = 2 
control Ch 2 < 
potentiometer 
> 
Fig. 1. 
it- 
i 
Fig. 2. Fig. 3. 
5. The resistor Ry and capacitor Cy are fixed. 
The equations corresponding to this model (for U = 0) are 
(1) 
= k,AU,, 
di 
L ++ Ri, = kU 
+ Ri,= 0 
di 
L' + R'i,= 
for U, < 0, (3) 
4 Rig=— | 
dU, 
GP + Uy, =U (ir, in), (4) 
AU, = (5) 


Here AU; is the voltage on the grid of tube 1, L and R are the inductance and resistance in the plate circuit 
when tube 2 is working, L* and R’ are the same when tube 2 is not conducting, ky and ly are the gains of the first 
and second stages ;k, = Uymax/AUymax, ke = F*¢Uymax/Uymax; b is the feedback coefficient, Up, is the feed- 
back signal and Us, is the voltage supplying the feedback circuit. Function @(i,,i,) can be put in the form 
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whic 


1 for M(i,)=1, D(i,) = 0, 
i:)={ for -0, © 
—1 for ®(i,)=0, O(i,) = 


Here &(i,, i,) has been put in this form on the as- 
sumption that Py and P, do not operate together. 


We introduce relative quantities 
max’ Us max’ max 


derive (0), (03), og) and eliminate ¢ and v; 
the system of equations then becomes 


’ 
7, + =0 
do, 
6, A +4,=0 
iq 
Fig. 5 1 for D(2,)= 1, = 0, 
M(o,, 5) =, 0 for D(s,)=0, 0, (10) 
for M(s,) = 0, D(o,) = —1. 
Figure 4 shows (04, G3): 
1 for ® = 1, 
(11) 
f for < (249) = 0, 
for (cy) = —1, 
={ O for (12) 
Here 0 yp, Oge are initial values, and 
= = = is L = = = 


This set of equations has been put down with due regard to the nonlinear functions Fy(y) and F,y(y%), 
which have the forms 


F, (19) = 0 for »>0, (13) 
for 
= \ 0 for <0. (14) 
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Here F* (U;) 
F, (49) = — =— 
F 1 (U1)max 
(U1) 
F, (49) = — F,(v) = 
2 (v) F\(Uvmax’ 
FY 1)max = 41 maxi, 
= 
whe 
Hence, the loop is specifiable in terms of the six parameters T,, Ty, Ts, y, €1, €- 
Figure 5 shows the structural diagram of a dynamic analog to the loop. There are three aperiodic links, on 
and four nonlinear characteristics. Phase-space and point-transform methods are best used with a nonlinear sys- P j 
tem of this type. The units in the loop have piecewise-linear response, so the phase space splits up into several ic 
regions, in each of which we have a set of linear differential equations, whose right halves differ from one region 
to another. The equations for the boundaries are known, and so we can follow the phase loci throughout the 
space, by transforming one surface to another point-by-point. We can establish the structural bifurcation in the tane 
phase space exactly by point-transform methods, for different ratios between the parameters, and can thus find viat 
what types of oscillation can exist. 
The system (7)-(12) shows that the phase space is three-dimensional (04,03, n). The equations for n > 0 = 
and n < 0 differ, so the plane n = 0 cuts the space into two parts symmetrical about the origin. The 04,02, 9 
phase space is bounded by the planes o, = 0, o, = 0, n = 1, n = —1, og = Y, Og = y, Which minimum and max- as pl 
imum values may be taken as being variable. The phase space is divided into several different regions by switch- ferer 
ing surfaces, because ©(0,,0,) is nonlinear. These surfaces are planes, and so the space is cut up into several and ; 
rectangular parallelepipeds, which overlap because ®(o4, 0) has regions in which it is not single-valued. Figure and « 
6 shows the phase space. The internal region is bounded by the planes oy = €, and og = é, because from &(o4, 43) O,= 


we have excluded the region in which 4(o,) and #(o,) are, for all values of 0, and og, simultaneously equal to 
+land —1, 
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Consider the region n > 0 in the 04,03, space. There are five regions where the functions are not single- 
valued, namely derhjxyl, where (04,09) can take the values —1 and 0; fribvulw, where (04,03) can take the 
values +1 and 0; and chrfuzjl, where ©(o 4,02) may take the values +1, 0 and —1. There are also regions 
cfm'mqzuq’ and chk'ktzjt’ where (04,03) can have the values +1 and —1, and so some initial deviations can 
give rise to cases where (a4) = +1 and (03) = — 1 simultaneously. 

Now (7)-(12) show that, when n > 0, og and n increase, whereas o, decreases; hence the phase loci can 
leave the volume hinmgjws via the planes iwjh and mqjh. In the latter case, the region where both relays are 
operated together must be reached. Hence we consider only the motion for initial deviations such that the phase 
loci leave the yolume himnqjws via the plane iwjh. This region is the one to which applies the condition that 
both relays do not operate together. The bound to this region is the locus through hj. This limiting case will 
correspond to some definite reJation between the loop's parameters. We may delineate a region of parameter 
values for which the two relays are not operated together, for a given set of initial deviations. 


We eliminate time from (7)-(12) and integrate to derive the surfaces in phase space that intersect to give 
the phase loci: 


when n < 0 
4 — (ce) 
(cs), 
— (a) 
when n > 0 
1 (2)| + (17) 
[y — + yD(o), 
(18) 


“= 
where a = b = Ts/T , O49, Ogg, no are initial values and = (04,03). 


It is clear that the loci are found as the intersections between cylindrical surfaces whose generators are 
parallel to the oy and og axes. The point O is a stable special point. Then the projection of a locus on the o,,y 
plane never depends on oy, while o, defines in the o3,n plane only the point where ©(04,03) intersects it. The 


projection on the o4n plane never depends on og, while oy defines in the o4,y plane only the point where 
(04,0) intersects it. 


Hence, our problem becomes to study the motions of the projections on the ¢4,n and dg, planes simul- 


taneously, because these projections represent fully all possible types of motion caused by the various initial de- 
viations in the 0,03, space. 


The Phase Planes for the Loop 


We consider the 04, and og, in the o4,09,n space as the phase planes for the loop. We represent them 
as planes of many sheets, because it is convenient to project each of our rectangular parallelepipeds on a dif- 
ferent sheet. The projections of the loci on the o4,n and o 3,7 planes are given by (15)-(18). Each of the oy,» 
and 0», planes consists of six sheets. The lines of intersection are o4 = €4, 04 = €g, n = 0 on the oy,y plane, 
and Og = €4, Og = €:, n = 0 on the Gg,n plane. The oy,n and og,n planes have bounds a, = 0, og = 0, n = +1, 
O1= 7. 
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Figure 7 shows the 04,n and o2,n phase planes, and indicates how the sheets are denoted. Each plane has 
two projections of the special point O, which is the equilibrium point for the system. The other special points 


' Me outside the phase planes, on account of the restrictions that have been applied. Some of the volumes in 


0 Og, Space overlap (Fig. 6), so sheets I, III, IV and VI in the o4,n plane, and I’, III’, IV’ and VI" in the og,n 
plane, have identical bounds, and are exactly superimposed. 


Fig. 7. 


The condition that the two relays do not operate together limits the parameter space available. The re- 
gion of initial deviations for which this condition must be complied with must be given when the bounds are to 
be found. Here we take this region to be that bounded on sheets V and III" by the loci that pass through 0, = y, 
n = 0 and og = €y, n = 0. The equations for these loci give us the relation between the parameters that defines 


the bounds in the following form: 


We may analyze the 04, and o2,n phase planes to find two main regions of initial deviations; to these 
will correspond different types of motion. In the first, the tendency is always to approach the point of rest. 
Within this region the image point always approaches the rest point; the transient may or may not show an over- 
shoot. There will be two types of motion within this region. 


1, In the first type the point lies either in n < 0 or in n > 0, and reaches O without passing outside its in- 
itial half-plane. On the a, or o, axis there is a bound R to the initial deviations o, for which such motion can 
occur, 

Figure 7 shows that this motion will occur in the top half-plane if the change-over point A‘ on o2,n lies 
within the region bounded by the locus M* that touches the line a, = €g. To find R on the oy axis it is, however, 


better to consider the points where M' and sheet I’ (which passes through A’) intersect the n axis (n,, and n°). 
The bound R corresponds to n,,, = n*. We use the equations for the loci on sheets V, III’ and I’ to find an ex- 


pression for R: 
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R is found on the og axis by putting o yg instead of ogs in (20). 


(20) 


2. In the second type the point lies at first either in n > 0 or in n < 0, and passes to the other half-plane 
before approaching O. On the oy or oy axis there is a bound S to the initial deviations for which such motion 
exists. Figure 7 shows that, if such motion starts in the top half-plane, we must have the point A in the 4,7 
plane to lie in a region bounded by the locus M, which touches the line a; = €y. Bound S on oy will correspond 
to Ron og. We use the equations for the loci on the upper half-plane, with (20), to get this expression for S: 


| | 


f is the ordinate of the point F. 


The bound S divides off the region within which the rest point is approached. This second type of motion 
may not exist for some ratios between the parameters. Then S vanishes, and the bound R takes over its function. 


The second region is the one in which the transient is oscillatory. The motion of a point initially in y > 0 
or n < 0 is then such as to take the point in turn through all six sheets, and back to the initial sheet, The loci 
may spiral inward, so that the rest point is eventually reached, but closed loops may exist. We shall smdy the 
motion in this region by point-transform methods. 


The Point Transform 


To seek limiting cycles and to study loci in the o4,n and og,n planes we examine the point transforms of 
a straight line into a straight line on each of the planes. It is best to see how the points on the o, axis transform 
to points on the same axis in the o,n plane. The point transform for any point is made up of two correspondence 
functions. For o, these function: are A, which defines the transform from points on the oy axis with o; > €,4 to 
ones with €; > o,> 0 (top half-plane), and B, which defines the transform from points on the oy axis with ¢,; > 
> 04> 0 to ones with o, > ¢€,. Function C defines the transform from points on the o, axis with €,; > a, > 0 to 
ones with o, > €;, while D defines the transform from points on the o, axis with o, > €, to ones with ¢, > o, > 0. 
The two phase planes are the same, so A is the same as D, and B is the same as C. 


The point-transforms diagram will consist of two parts, each of which defines the transform for the coor- 
dinates that correspond to itself. The final and initial o, and og values are related, so the diagram will be a 
space, and the functions will define surfaces. We must consider both parts in order to follow the motion of the 
image points as the number of transfers increases. The fixed points can be found by considering one coordinate 
only, because the functions are the same for both parts. We consider the point transform for o, below, and ex- 
amine the behavior of A and B in general terms. 


Half-plane y > 0. Let the coordinates at the start be o, =u, n = 0, og = w. We consider A, so u> €,, 
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id O<w< Now take some time t® = ty + ty + ty, where ty, 
t, and ts are the times spent by the image point in sheets V, 
YUNA IV and VI, when the coordinates are o; = v, n = 0, Og =F. 


V4 ; We use the solutions to (8) and (9) for ty, t, and ts, and 
7 eliminate from these ty, t,, tg, w, and intermediate values of 
A certain coordinates; we then have relations between u, y and 3 


> b 


V, r=(%+7~ 


7 ™ Half-plane  < 0. Let the coordinates at the start be 

Fig. 8. 04=V, Og =r, n = 0. Consider some time t** = ty + ts + tg, 

where ty, ts and t, are the times spent by the image point in 

sheets III, I and Il, when the coordinates are o, = Uy, Og = Wy, 

n = 0. We use the solutions to (7) and (9) for ty, ts and ts, and eliminate wy, ty, t; and tg, and certain intermed- 
iate values of coordinates; we get expressions relating uy, y and r: 


(24) 


where ¢ is the ordinate of the point where the trace of sheet I meets the line o, = €g. The Appendix examines, 
in a general way, the behavior of A [which is defined by (23) and (24)] and of B [which is defined by (25) and 
(26)] for the case a > 1, b> a. 


Figure 8 shows part of the point-transforms diagram for oy in u, v, r coordinates. Here A (a meet is 
shown for 0*v/dr* < 0, and B for 8*v/dr* > 0. The second part of the diagram relates to og, and is in u, w, r 
coordinates. The A and B surfaces meet in one line LL’ only. Figure 8 also shows that, with this deposition of 
the surfaces, u will change as the number of transfers increases in such a way as to approach LL' (whatever r), 
and hence that the fixed points defining the limits to the cycle must lie on LL'. Now limits can exist on both 
planes only if u = r, so the point corresponding to the limit will be the point where LL’ meets the plane u = r. 
Consider the plane v = const in order to see how LL’ lies, with respect to the plane u=r. This latter plane will 
give a trace of the form u = r on the r = const plane, whereas LL’ will give a projection Ul’, whose equation is 


given by (23)-(26) as 
an 


Now the equation for ll’ is complex, so we cannot consider its behavior fully, but some general deductions 
can be made. The derivative du/dr for Jl’ is given by 
a 
=) 
du a(yf—e)(u> — 
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dr wher 


1 1 a-1 


Now du/dr is positive for all r, and when r = rg we have du/dr = oo. The expression for du/dr* is very 
complex, and even the sign of the derivative cannot be found from it. We may say, though, that if r is close to 
Io, We must have @u/dr* negative, because du/dr, while positive can be less only than +00. 


The phase~plane diagrams also show that the end of Jl” (i.e., 2") must always lie above the straight line 
u =r, since the bound to the point transform, r¢, corresponds to f = 1, and hence, from the transform for A, to 
u=oo. Butr =r* gives us for the transform of B a value of uy less than re. The shapes of the A and B surfaces 
show that ll’ may, in general, have several inflections. 
Bifurcation Relations for the Parameters 

The following disposition of Ul" relative to the line u = r can occur. 


1, Line I" lies above u = r (Fig. 9). The motion is then stable for all initial deviations, 


2. Line UW" cuts u = r, but does not reach the bound r = r, (Fig. 10), There can then be any even number 
of intersections, to which will correspond stable and unstable limiting cycles. Stable motion will still occur for 
some initial deviations. The dividing line between the first and second cases corresponds to UJ” touching u = r, 


Then we must have u =r and du/dr = 1, We use these conditions, with (27)-(30), to give the bifurcation relation- 
ships as 


i—a 
Y_ __ (31) 


i—a 


! 
' 
% 


> 


Fig. 9. Fig, 10, Fig. 11. 


Now both conditions, u = r and du/dr = 1, are allowed for in (32), whereas only u = r is allowed for in (31), 
so the roots of (32) are also roots of (31). The two equations give the critical relation between the parameters 
when they have a common root. If (31) has extra roots, Ul’ and u = r will cut as well as be tangential. If (32) 
loses or gains a root, or (31) two, we get a change in the type of motion, and may find a bifurcation relation. 
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v Y ZY 3. The line 22" cuts u = r, and J passes through the entire re- 
Y GY gion r, =r <r* (Fig. 11). There :nust then be at least one inter- 
Z c 44). 
is Stabi |Z section of with u =r. The A and B surfaces considered here are 
GY regi so placed that this point corresponds to a stable limit. All devia- 
Z all tions within the second region will then give unstable motion. 
w—7Z The boundary between the second and third cases occurs when 
%Y l lies at the point u = rg. This point will correspond to cycles that 
Y include parts of loci M and M’. The coordinates and equations for 
“7 the loci give us the bifurcation relation for the second and third 
types in the form 
where 


G= (a — 1) 


Hence our study of the possible positions of il’ relative to the line u = r on the v = const plane has given 
us equations from which we may delineate areas in phase space with differing structures within the region where 
the transients are oscillatory. There are three main regions in phase space. In the first, the motion is stable, in 
the second, the motion may or may not be stable, and in the third, we get hunting. The boundary between the 
first and second regions is the stability boundary, and is defined by (31) and (32). Equation (34) defines the 
bounds to the third region. 


The loop must always be stable, and so in any real case we have to find only the stability limit, i.e., to 
find only that relation between the parameters for which (31) and (32) have a common root. 


_ Figure 12 shows the stability bounds in the a,b plane for a loop with y = 1, €y = 0.1, eg = 0.2. The projec- 
tion of 21° for a from 0.1 to 0.9 is of a shape such that only one contact between it and the line u = r is possible 
in the vy = const plane; hence (31)-(33) give us the relation that defines the stability limit. 

APPENDIX 


We examine the shapes of the surfaces defined by (23)-(26). We assume that a < 1, b > a, which is so for 
the loop in an electrical autopilot. 


Case a < 1, b> a. We consider the sections of the surfaces formed by planes having r, u and u, constant. 
The surface defined by function A has derivatives dv/du and 07v/du* of the forms 


b+1 
fre,” (35) 
(1+ f)?(ul?— 
+ 1) 
b(1 + (ul? — el? 


(36) 


Now (23) and (24) show that, when u varies from ¢, to o, v varies from ¢, to v* = €;(f/1 + f)>. Then 
is negative for all u > whereas 0*u/du" is positive. Also and 0*v/dr* have the forms 


(4 + fy (37) 
a (yf — ea) (ul — el”)? 

(ul? — (yf — 
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wh 


These expressions show that, if r > rg, 1.e., if f > €,/y, Ov/dr is positive, The latter is +o ifr =r. Then 
Tc is the bound to the point transform, and corresponds to fo = ¢/y. Now dv/ér" is negative when b <1 for all r. 
If b > 1, it is possible for &*v/ar* to change sign once or twice, since [b — 1 — f(1 —a)] (yf — &) — yf” — yf can 
become zero for two values of f. 


Our study of surface A thus shows that this surface intersects a plane r = const to give a falling curve for 
all r such that rc < r < r®, and that y will tend to a definite limit if u-»+ a. The sections in a plane u = const 
depend on b, and can be of three types. 


Let us study surface B (function B). The expressions for dv/du, and y/auj are 


a 

ou? ag’ (yp — (40) 


Now function B has a bound uy = Ugg, which corresponds to ¢, = ~ €/y. 
Then (25) and (26) give «'s that u; changes from Ujc to u% = r as g varies from ye to ~1, while y varies 


But dv/du, is positive if gy < go. However, #v/duj is negative when y < gp. 


Let us examine the section of surface B by planes with u, = const. The expressions for dy/@r and dv/ér* 
can be derived from (26), but they are very complex, and the signs of the derivatives cannot be found. We use 


the equations for the loci in the lower half-plane and the solutions to (7) and (9) to relate y and r in an expanded 
form: 


ou = — gel, = gels, 


where c = 1/Ty, h = 1/T3, and ng and oy are the coordinates at the end of the fourth (t,) part of the motion. 
Then dv/dr can be put as 


Or Or Borg Oty Ong 
The partial derivatives in this expression are given by 
i dv 


i 
cect fey ofa — Ft = 


These formulas show that 0n4/dr is negative, that dv/dn, and dv/d0 y are positive when < ~@/y 


that d0 y/dt, is negative when y < —€,/y and ts > 0, and that dt,/dn, is also negative. Hence dv/ér, which 
can be put as 


ay 
av — a Y 
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from v* to €y, where 
— i, 


is negative for all the r and u within the set range. 


Now @v/ar? is given by 
x[= uy’ (1+ o—a) (43) 


The expression in square brackets in (43) differs from the similar one in (42) in that the negative term is 
multiplied by a quantity greater than unity which decreases as r increases. Hence, #v/ar® can be either positive 
or negative, and the change between the two may occur once. The reason is that dv/dr is negative for all r. 


Our study of surface B shows that the sections by planes with r = const are rising curves, and that when 
Uy + @, r-+ 00. The sections given by planes with u, = const can take three forms; which form is found depends 
on u, r and the relation between the parameters. 


I wish to thank N.A. Fufaev, who made some very valuable suggestions about the content of this paper. 
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COMPENSATION OF A CONTINUOUS AUTOMATIC CONTROL SYSTEM 
BY MEANS OF A DELAY ELEMENT FILTER 


Wang Sin-Min 
(Moscow) 


Methods are considered for the design of a delay-element filter for compensating 


a continuous automatic control system. The filter parameters are determined for a few 
simple cases. 


For some time past papers have appeared in the literature, dealing with the use of a delay-element filter 
(DEF) for improving the dynamic quality of continuous open-loop and closed-loop systems, When compared to 
an ordinary compensating network a DEF has several advantages: flexibility in the selection of its parameters 
and in adjustment, convenience of synthesizing the filter from a given pulse characteristic, etc. These advantages 
make filters of this type particularly convenient for the compensa - 
tion of systems and for use in simulation [1-6]. The characteristics 
of the DEF also permit one to design a system which attains a 
steady-state mode in a finite time interval (1, 2, 3, 7] and to com-+ 
pensate the delay in an automatic control system [8]. 


In this paper, the use of a DEF as a compensating device in 
continuous automatic control systems is considered. 


The DEF shown in block-diagram form in Fig. 1 has a trans- 
fer function of the form 


Fig. 1. 


D(p) = + aye~ ase — 2PT +... age PT, 
where T is the delay time of each delay element and n is the order of the filter, characterizing the number of 
delay elements. 


The problem of compensating a continuous automatic control system by means of a DEF consists of the de- 
termination of the order of the filter and the values of its parameters for which the system takes on the required 
properties. It is preferable to choose the order of the filter as low as possible, and, at the same time, its param- 
eters must be optimum in a particular sense. 


Determination of the Order of a DEF 


Let us examine the automatic control system, the schematic of which is shown in Fig. 2, where the expres- 
sion 


represents the transfer function of the permanem part of the system, while D(p) represents the transfer function 
of the compensating DEF. A considerable number of networks corresponds to the form of transfer function given 
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M, M, by K(p). Moreover, the compensation of a system by such a circuit 
Dip) M(p) represents a more difficult case then the compensation of a system by 
a circuit having a polynomial in p in the numerator of the function 
K(p). 


Fig. 2. Ideal compensation of a system requires that for a static system 
the relationship 


D(p)K(p) = C (1) 


be satisfied, while for an astatic system the relationship 


D(p) K(p) = C/p (2) 


be satisfied, where C is a constant which should be chosen as large as possible. Because of the limitation of per- 
missible power and the existence of saturation in practical systems, condition (1) or (2) can be satisfied over a 
definite frequency range only. 


We know from the papers which dealt with the case of the compensation of an open-loop system [1,2] that 
condition (1) can be satisfied over the range of frequencies of interest to us by making use of a DEF. In this case 
the sufficient order of the filter is equal to the order of the system being compensated. 


It follows from what has been said that to compensate a static, automatic control system one may select 
the order of the filter equal to the order of the open-loop system, while for an astatic system the order of the 
filter will be one lower. 


In a real system, the time constants of the fixed part of the system do not all have the same effect on the 
properties of the system. It is obvious that to ensure the required properties of the system it suffices to take into 
account those time constants which have the greatest effect on the system's properties. Hence, the order of filter 
can be chosen tentatively to be equal to the number of these time constants. 


Having chosen the order of the filter, we will examine the question of choosing its optimum parameters. 


Optimum Parameters of a DEF 


The problem of obtaining optimum parameters is equivalent to the basic problem of the theory of the ap- 
proximation of functions in mathematics. Different criteria exist for the selection of optimum parameters of 
controlled systems. We will examine here the method of selecting parameters, which is based on the theory of 
the approximation of a function by the method of the expansion of the frequency characteristic of a control sys- 
tem into a Taylor series. The parameters are called optimum when the frequency characteristic of the compen- 
sated system deviates by the least amount from its value when w = 0. 


The method of expanding a frequency characteristic into a series was first used by Braude in the synthesis 
of reproducing systems [9]. The use of this method for the choice of optimum parameters of control systems was 
described, for example, in papers [10,11]. Below we will examine the application of this method of choosing 
parameters to our case. 


The transfer function of the closed-loop control system shown in Fig. 2 is equal to 


K (p) D(r) M (p) 3 


M (p) = + 
A (p) = by + ++ +++ bmp™ + dy + aye—T 


If we expand the square of the modulus of the frequency characteristic | Ke (jw)|* = 4 (w) into a Taylor 
series about the point w = 0 


(@) = +9" (0) 4... (4) 


(the odd derivatives are identically equal to zero) we obtain the conditions for optima [10): 
(0) = x (0) = = 0, (5) 


In this case the frequency characteristic A (w) has its smallest deviation from A(0). The problem of choos- 
ing optimum parameters consists of finding expressions for the conditions for optima (5) from the filter parameters, 
appearing in Eq. (3), and then determining their values. 


We find ar expression for A (w) from Eq. (3): 


Expanding the numerator M (jw) M(~ jw) into a series (the odd terms are identically equal to zero) we obtain 


M (jw) M (— jw) = My + Myw* + Mw +---, 


al 
2 1) (/) d (jw)* d (—ja)**—* (1=0, 1, 2, ds 


Since here the differentiation does not depend on the sign of w then, after replacing jw and —jw by p, we 
obtain 


al 
Mu= X(—1) | (8) 
Analogously, for the denominator of expression (4) we have 


A (jo) A (— jo) = By + + 


al 
kt (p) A@!—*) (p) 


The expression for 4 (w) can now be written in the following form: 


M, + + Myo +--- 


It should be noted that relationships (6) and (7) permit one to calculate Mgz and Ag; directly from the 
transfer function K,(p) without having to find the square of the modulus of A (w). 
Let us examine the expression for the deviation of the frequency characteristic A (w) from (0): 


— Mods) w* + (Mid, — Mp Ay) 


where 
where | 
4 
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A comparison of expressions (4) and (9) shows that the quantity A"(0) will be equal to zero if the condition 
M,/Ag Mg/Ag is satisfied. 


Analogously, the derivative 4“) (0) will be equal to zero if the conditions Mg/A, = Mz/Az = Mo/Ag are 
satisfied, etc. Hence, the conditions for optimum filter parameters can be written in the following manner: 


My 10 


In the case of an astatic system (Mg = Ag) expression (10) takes on the form: 
My = Ay (2 = 1,2,...), (11) 


In determining the optimum parameters (a9,a4, . .. ,an, T) of a DEF from Eqs. (10) or (11) one must use the 
n + 2 equations which correspond to the n + 2 first (even) derivatives of the square of A (w) at the point w = 0 


being equal to zero. 


In this way, the selection of the optimum parameters of a DEF consists of the following: values of Mg; and 
Agi are determined from expressions (6) and (7) under the conditions 1 = 0, 1,...,n,n+1,n+2, where n is 
the order of the chosen filter; all the ratios Mjz /A2; are equated to the value of Mg/Ag (for a static system) or 
to unity (for an astatic system); all the desired parameters of the filter are determined from these equations (con- 
ditions for optima). 


In determining the filter parameters one can take into account other restrictions, such as the required values 
of the error coefficients, limiting values of several filter coefficients a4, etc. 


The transient characteristic of the compensated system must be checked with the chosen filter parameters. 
If the characteristic does not satisfy the required properties of the system, then it is necessary to raise the order 
of the filter and repeat the calculations. 


Comparison of Ordinary Filters with a DEF 


The compensation of a control system by a DEF is achieved by introducing, beside a signal proportional to 
the error at a given instant of time, signals proportional to the errors at previous instants. This, in actuality, is 
equivalent to introducing into the law of compensation, derivatives of the error. But, in comparison with conven- 
tional differentiating devices, the described filter is more stable in the presence of interference. 


Thus, let us compare three differentiating devices from 
db this point of view: an ideal differentiating circuit having the 
transfer function K,(p) = Typ + 1, a passive differentiating cir- 


Vall _ __Typ+1 
; cuit with transfer function K,(p) = 0.05Typ + 1 and a DEF with 


an transfer function D(p) = 10.5-9.5e79-1T yp. Over a certain range 
of frequencies the characteristics of all three devices practic- 
ally coincide. 


S70 200 307000 20 Figure 3 shows logarithmic plots of their amplitude char- 
acteristics. 


It can be seen from Fig. 3 that in the case when noise is 
a component of the input signal, the DEF (curve 3) will ensure 
a better result than will the ideal differentiating circuit (curve 1) and the passive differentiating circuit (curve 2). 
If the noise spectrum is concentrated in the vicinity of the frequency wp, then the delay time of the delay ele- 
ment T should be chosen by means of the following relationship: wnT = 2kn, where k is a positive whole number. 
If T is selected in this way the noise passes through the DEF without the increase which occurs in ordinary differ- 
entiating devices. 


The DEF differs from other compensating networks in its flexibility with respect to changes in parameters. 
Calculated optimum parameters of a compensating unit do not always yield an optimum control process in a 


Fig. 3. 
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practical system, since an idealized transfer function of the system is used in the calculations. Consequently, 
the filter parameters are made more precise by experimental means during alignment. This is yery simple to 


achieve in a DEF by changing the settings of certain potentiometers. 
Examples of Compensation of Automatic Control Systems 


Example 1. Given a system, the circuit of which is shown in Fig. 4, it is required to determine a compen- 
sating DEF with optimum parameters, 


We choose a first-order DEF with transfer function 


D(p) = “pT, 
(p) = ag + age 
y The transfer function of the open-loop system will be 
Fig. 4. W (p) = 
Let 


T= = Ag, Ky Ti = A. 
Substituting these relationships into the expression for W(p) we obtain 


Ay + 
CES) 


The transfer function of the closed-loop system will be 


A, + Aye a A(q 
K = = » 


Let us find the derivatives of the numerator and denominator of the transfer function K,(q) when q = 0. 


Derivatives of ; Derivatives of Derivatives of Derivatives of 
numerator — denominator numerator denominator 
d0A aA | 
= Ar age = 40 +41 age ag = — 
= AyT* ag =2+ A;T* 


From relationships (6) and (7) for finding values of Mgy and Ag; and from the conditions for optima (11) 
we obtain 


1 + T) = 0 when = 2, 
whenl =4. 


Solving this system of equations, we obtain | 


3 


4 

i 

q 

4 

| 
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0 Bw 0 sec 
Fig. 7. Fig. 8. 
The amplification factor is 
K = (a9 + = 1/Ty(Ag + Ay) 
or 
| 3 


Variations of Ag, Ay and KT, as functions of the relative delay time of the filter T are shown in Figs. 5 and 6. 


Starting with the required value of the amplification factor K, we find the delay time of each element of 
the DEF from Fig. 6. The corresponding values of Ag and A, of the filter can be determined from Fig. 5. 


For example, if the controlled object has Ky = 1 and T, = 1, while the required amplification factor K = 3, 
then from Figs, 5 and 6 we find T = 0.36, ay = 9.9, a = — 6.9. 


The transfer function of the chosen filter will be 


D(p) = 


426 


| 
4 
Fig. 5. Fig. 6, 
\b 
ol 


= % The amplitude-frequency characteristics of a compensated 
os Dip) , system (curve a) and an uncompensated system with the same 


pli amplification factor K = 3 (curve b) are shown in Fig. 7, The 
transient characteristic of the compensated system turns out to 
Fig. 9. be slightly oscillatory (Fig. 8). 

Example 2. Given a system, the circuit of which is shown 
| in Fig. 9, 


optimum parameters. 
We select a second-order filter with transfer function 


D(p) = ap + aye "PT + age™*PT, 
The transfer function of the open-loop system will be 


(ao + + age—*P7) 
p(T®p* + + 1) 


Let pTy = 4, T = TyT, KyTyq = Ag, KyT = Ag, KiT gay = Ap. 
Substituting these relationships into the expression for W(p), we obtain 


W (p=) 


Ao + +. 


= ~~ + +1) 


The transfer function of the closed-loop system will be written in the form 


Ay + 4 


K, 
+ +9 + Ag+ + BCA)” 


(9) 


We find the derivatives of the numerator and denominator of the transfer function K,(q) when q = 0. 


| 
Ay + Ai + As Aa + Ant 
= + =1 —(A;—2A,)T 
(A; + 16A,)T* = (A; + 164,)T* | 
+ 64A,)T* = (Ai + 


From relationships (6) and (7) for determining Mgz and Agj, and from the conditions for optima (11) we i] 
obtain 
+ (46 + 27) Ay + (46 +47) for 2, 


(67 + 6ET? + T*) Ay + (127 87%) Ay 6 — 1289 for = 4, 
(207? + 4+ Ay + (1607°4160ET* + 327") 60 for 6. 


1,4, 
| 
\ 4 
= Fig. 10, Fig. 11. Fig. 12. 
az a4 Qs ‘tf. a7 aé ag 10 
Fig. 13. 
la 
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Fig. 14, 


The amplification factor is 


+ a4 + ag) Ky = + Ay + Ag) 


KT, = Ag + Ay + Ag. 


The variations of Ag, Ay, Ag, and KT, as functions of the relative delay time T of each element of the 
filter, for different values of €, are shown in Figs. 10-13. Starting with the requized value of the emplification 
factor K one can, making use of these curves, determine all the filter parameters, 


For example, if for the controlled object Ky = 1, Ty = 1, € = 0.4 and the required amplification factor K = 1, 
then making use of Figs. 10-13, we obtain T = 0.5, Ag = 4.71, ay = ~ 6.90, ag = 3.17. 


The filter's transfer function will be 


D (p) = 4.71-6.90e*"P + 3,17e7P. 


The amplitude-frequency characteristic of the compensated system is shown in Fig. 14, while the transient 
characteristic is shown in Fig. 15. An uncompensated system with an amplification factor K = 1 is unstable. 


SUMMARY 


1, Delay element filters make it possible to improve the dynamic properties of automatic control systems, 


2. Optimum filter parameters may be chosen by means of the method of approximation by an expansion 
of the frequency characteristic into a series (Braude's method). 


Other restrictions imposed on the system, such as the specification of particular values of the error coeffi- | 
cients, can be taken into account when determining the filter parameters, i 


3. For the case of the compensation of simple closed-loop systems, including circuits with transfer functions 


Ks 
+1) p (Tip? + ip +1) 


formulas and curves have been obtained for determining filter parameters. 
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INVESTIGATION OF THE CHARACTERISTICS OF PNEUMATIC JET ELEMENTS 


L.A. Zalmanzon and A.I. Semikova 


(Moscow) 


The article deals with an investigation of pneumatic jet elements, the operation 
of which is based on the utilization of the characteristics of pressure changes obtained 
when the impact tube is moved through different sections of the jet flowing from a 
nozzle. Examples are given of the possible applications of such elements in pneumatic 
devices. 


Pneumatic devices are built at present, as a rule, with nozzle-vane type elements [1]. In the solution of 
many problems which arise in the field of pneumatic devices, the application of jet elements, henceforth referred 
to as nozzle-pipe elements, can be more effective. Their operation is based on the utilization of the character- 
istics of the distribution of pressure in an air jet. 


Nozzle-pipe elements can be divided into two types. Elements in which the Pitot tube moves through the 
section of the jet emerging from the nozzle belong to the first type (Fig. la). To the second type belong ele~- 
ments in which there are no moving parts and variations in 
pressure on the pipe are brought about by the deflection of the 
air jet emerging from one of the nozzles, by a jet emerging 
from the other nozzle (Fig. 1b).* 


\) eon Sy \) Jet elements have several advantages over elements of 
— the nozzle-vane type in that there are no precise, ground de- 
WR tails in them, it is not necessary to have hermetically sealed 
2 connections in them, and they work more reliably over a wide 
eR temperature range. A disadvantage of jet elements is their 
ay ne relatively high consumption of air. This last fact is insignificant 
if powerful sources of air exist in the controlled object itself 
(gas turbine and uniflow air-compressor motors, ergazlifts in the 
petroleum industry, etc.). In pneumatic equipment designed 
for general industrial application, the consumption of air can be reduced at the expense of a reduction in work- 
ing pressures. 


The use of nozzle-pipe elements of the first type is expedient in those cases when it is required to match 
variations in pressure with mechanical displacements. Nozzle-pipe elements of the second type are purely aero- 
dynamic and serve to vary one pressure as a function of other pressures, In this paper nozzle-pipe elements of 
the first kind are investigated. In part 1 examples are given of their possible application in pneumatic devices. 
The following parts are devoted to the theory of these elements. 


Fig. 1. 


The characteristics of variations in pressure occurring on the receiving tube as it moves through different 
sections of tie air jet emerging from the nozzle are investigated. Specific problems arise in this section, in con- 
nection with the small dimensions of nozzles in nozzle-pipe elements, and in connection with the possible effect 


*Concerning the use of jet elements in pneumatic devices see also [2]. 
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of the tube dimensions on the obtained characteristics. Section 2 clarifies the extent to which one can, in cal- 
culating the characteristics of nozzle-pipe elements, make use of data related to free jets. The analysis of 
functional transformations carried out by elements of the type under consideration is given in part 3. In section 
4, possible divergences between the characteristics of nozzle-pipe elements and the characteristics of free jets 
are examined further. Here consideration is given to questions having special significance in cases where the 
relative dimensions of the receiving end of the pipe are so large (in connection with the small dimensions of the 
nozzle opening) that they cannot be considered negligibly small in comparison with the dimensions of the jet's 
cross section. 


1. Examples of the Possible Use of Nozzle-Pipe Type Elements in Pneumatic Devices 


Nozzle-pipe elements of the type being examined can be used effectively in devices for the transformation 
of mechanical displacements into air-pressure variations and in devices for the reverse transformation. 


In the transformation from mechanical displacements into changes in air pressure, the problem of obtain- 
ing linear characteristics is sometimes set; more often, however, the corresponding device must achieve functions 
of a nonlinear converter. 
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The simplest converter is a single nozzle-pipe element. By displacing the pipe in a definite manner 
through different sections of the jet one can obtain different shapes of nonlinear characteristics (Fig. 2a,b,c). 
The shape of the characteristic P = f(s) or the characteristic P = f(h) depends on the size of the nozzle and pipe, 
on the input (supply) pressure, and on the initial placement of the pipe in a section of the jet. For example, by 
varying the input pressure and the placement distance h, one can obtain the family of curves P = f(s) shown in 
Fig. 2d. An element which yields such characteristics can serve, for example, as a model of the controlled ob- 
ject in the simulation of optimalizing control processes.* By using near to linear portions of the characteristics 
indicated above (see, for example, regions ab in Figs. 2a and b) one can carry out transformations which are 
practically linear, Using only one nozzle-pipe element one can carry out functional transformations for the 
achievement of which it is often necessary to use a complex to produce and finish cam mechanisms, lever mech- 
anisms, etc, 


. Nonlinear characteristics of a more complex form, for example, characteristics P = f(s) with two or more 
extremums (Fig. 2e) or others, can be reproduced by the combination of two or several elements of the type in- 
dicated earlier on. 


Nozzle-pipe elements can also be used to build a universal, nonlinear converter reproducing any given 


*In the simulation, the displacement of the maximum of the characteristic P = f(s) along the s axis can be pro- 
duced by changing the position of the supply nozzle, and along the P axis by changing the supply pressure Pjp. 
The static characteristics of the object are considered here. However, when a pneumatic chamber or a system 
of pneumatic chambers is connected to the output of the given element, the dynamic properties of the object 
can be reproduced on the model as well. 
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characteristic P = f(s). A characteristic P = f(s) of an arbitrary 
form can be obtained by the use of a set of nozzle-pipe elements 
with nozzles having an elongated cross section (Fig. 3a). 


For each of these nozzles the characteristic P = f(s) has a 
shape like the one shown in Fig. 3b for displacements of the pipe 
along the major axis of the nozzle’s cross section. If the move- 
ment of the pipe is restricted to the regions ab'c’ or b'c'd the 
resultant characteristics are those shown in Figs. 3c and d. As 
the input pressure is varied the level of the horizontal section 
changes and, accordingly, the inclination of the side branches of 
the characteristics indicated earlier (Figs. 3e and f) change. 


Let us examine the method of obtaining the characteristic 
P = f(s) having an arbitrary form, by means of a set of elements 
with nozzles of the above-mentioned form. Let, for example, 
the curve shown in the upper part of Fig. 4 be the specified one. 
Let us replace this curve by a piecewise-linear approximation. 
This we will construct in such a way that the projection of each 
piece of the approximation onto the s axis be equal to 5s (the 
value of 5s is determined from the characteristics shown in Fig. 
3e and f). The three middle graphs of Fig. 4 show elementary 
characteristics, each of which has a sloping region equal in length 
and angle of inclination to the corresponding section of the 
piecewise-linear approximation described above, Each of these 
elementary characteristics is obtained by means of one of the 
nozzles. Their required relative position along the s axis, shown 
in Fig. 4, is obtained by coordinating the position of the nozzles 
during their initial installation (Fig. 5) while the required inclin- 
ation of the sloping region of each of the elementary character- 
istics is obtained by setting a particular air pressure in the supply 
for the corresponding nozzle. The pressures received by the pipes 
of separate elements are added by means of a pneumatic, inter- 
throttle chamber which fulfills the function of a summator. The 
pressure at the output of the summator changes in accordance 
with a characteristic which corresponds to the given piecewise- 
linearized characteristic P = f(s) with, possibly, the only differ- 
ence being that the whole characteristic can be displaced by a 
constant amount of pressure, to one or the other side of the spe- 
cified one (see the characteristics shown by dotted lines on the 
lower graph of Fig. 4). The characteristic can be brought to the 
specified level either by the use of a nozzle-vane (Fig. 6a) or 
by means of an arrangement containing no moving parts and based 
on purely aerodynamic action (Fig. 6b). In the first case the dif- 
ference between the pressure supplied to the given element and 
the pressure at its output is created by the tightening of springs 
by a regulating screw. In the second case the functions of a dis- 
placement element is carried out by a pneumatic chamber analo- 
gous to the one used as a summator. As the pressure P;, varies, 
the quantity Poy follows the variations of Py, displaced by a con- 
stant value of pressure, determined by the adjustment pressuze 
Paqj- When the given adjustment element is connected to the 
summator, shown in Fig. 5, one obtains the characteristic P.+ = 


= f(s) displaced relative to the characteristic P}, = f(s) by a constant amount depending on Pg qj. 
When jet elements of the type being considered are used instead of nozzle-vane elements, the transforma- 
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tion of air pressure into mechanical displacement is also 
~~ more easily achieved. For comparison, Fig. 7a shows the 
Fig. 5 schematic of a converter in which a nozzle-vane ele- 


ment is used, while Fig. 7b shows the schematic of a con- 
verter using a nozzle-pipe element. In the latter case, 
nozzle-pipe 1 is used as a feedback element. As a result of the operation of the feedback loop there exists a 
strict correspondence between the value of the pressure at the input to chamber 2 (this same pressure is produced, 
as a result of the feedback loop, in chamber 3) and the variations in the distance between the nozzle and pipe. 
By locating the pipe in different sections of the jet one can design, even with relatively small nozzle dimensions, 
conversion programs for a relatively large mechanical displacement. The range of displacements can vary from | 
tenths of a millimeter to several millimeters or to values measured in tens of millimeters. | 


Fig. 8. : 


A study of the characteristics of jets shows that nozzle-pipe elements can be used effectively not only to 

carry out the above-mentioned functions, but also for other objectives. For example, making use of the fact that i 

for every given position of the pipe the ratio of the pressure on the pipe to the pressure in front of the nozzle does i 
not depend on the value of the latter (we have in mind excess pressures), one can use nozzle-pipe elements to 

carry out multiplication, division and squaring operations. Figure 8a shows a unit which illustrates the principle ; 

of the execution of these operations with nozzle-pipe elements.* When pressure P, is applied to chamber 1, the t 


*Compare this with the schematic of the computer presented in [3], which is used to carry out those same func- 
tions and which is made up of nozzle-vane elements. 
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same pressure P, is produced in chamber 3 as a result of the operation of the feedback system which includes 
nozzle-pipe element 2. To this feedback system P, is the supply pressure. On the same movable crosspiece on 
which the nozzle of element 2 is installed there is also the nozzle of another nozzle-pipe element 4, The di- 
ameters of the nozzles of both elements are the same as are the dimensions of the pipes. For every given value 
of h we have P,/P, = f(h) for element 2 and P3/P, = f(h) is maintained for element 4. Since, with the feedback 
loop operating, the steady~state value of h is the same for both elements, the values of f(h) in the right sides of 
the above-mentioned equalities are the same. Hence, P,/Py = P3/P,, whence it follows that P, = PyPs/P,. When 
P, = 1* we have Py = PyPs, when Pg = 1, correspondingly, Py = Py/P,. However, if P, = 1 and Py = Ps, then P, = Pj. 
The operations indicated above can also be carried out if one makes use of the characteristics of pressure varia~- 
tions which are obtained when the receiving end of the pipe is moved along the jet's cross section (Fig. 8b). If 
an extra feedback loop is introduced, the system shown in Fig. 8a can be converted in such a way that the opera~- 
tion of the extraction of a square root can be carried out. For example, for the circuit shown in Fig. 8c, where 
Pz = 1 and Ps = Py, we have P, = PyP3/P, = Pj or Py = VP. where P, is the pressure at the input and P, is the pres- 
sure at the output. 


Models of the universal, nonlinear converter and multiply -divide unit described above were built, using © 
nozzle-pipe elements,** at the laboratory of pneumohydraulic automatic control of the IAT AN SSSR. Experi- 
mentation with these arrangements confirmed the correctness of the design of their principal systems. 


The examined examples of the possible use of nozzle-pipe elements in pneumatic instruments can be ex- 
tended to hydraulic instruments. As it is shown below, the theory and design of nozzle-pipe elements are based 
on the use of the characteristics of free, turbulent jets, but these characteristics are common to air and drop li- 
quids. 


2. Investigation of the Possibility of Using the Characteristics of Free Jets in the 
Design of Nozzle-Pipe Elements 


The processes of the propagation of a free jet have been the subject of special investigations in aerody - 
namics (see [4,5]). Experimental investigations in these studies were basically carried out on the flow of air 
from nozzles, the diameters of which were measured in units of the ordex of s::veral centimeters or decimeters. 
Nozzles used in pneumatic devices can have diameters not exceeding tenths of millimeters in size, and only in 
certain special cases can they reach sizes of the order of several millimeters. In the investigation of the charac- 
teristics of nozzle-pipe elements in the role of primary elements in pneumatic devices there arises the question 
concerning the possibility of using the experience accumulated from the investigation of the flow of air from 
nozzles with relatively large aperture-cross sections. 


Let us first list some of the data concerning the free jets with the characteristics of which the characteristics 
of nozzle-pipe elements are subsequently compared. The profile of the longitudinal section of a free, axially - 
symmetrical, turbulent jet is shown in Fig. 9a. The initial region, of length hy;, contains a core of constant speeds 
(shaded in Fig, 9a) in which the axial speed of air particles is equal to their speed at the exit from the nozzle. 
The rest of the section in this region is occupied by a boundary layer. When h > hy in the main region of the 
jet, the boundary layer occupies all the cross section of the jet and the jet’s axial speed decrease; with increase 
inh, The decrease in the particle speed occurs for a constant amount of motion in the jet, and is dependent on 
the increase of the mass of the latter as a result of the drawing of air into the jet from the surrounding medium. 
The static pressure at all points in the free jet is assumed to be constant and equal to the pressure of the surround- 
ing medium. When flowing out into the atmosphere, this is the atmospheric pressure and, in that case, the pres- 
sure of the velocity head received by the pipe introduced into the jet appears together with the whole pressure 
measured in gauge units. 


The scheme for producing jets that is described above is a simplified one. In actual fact, the static pressure 
in the jet does not remain strictly constant. In the main region the intensity of the turbulent flow changes and, 
in accordance with this, one should divide the main region into separate zones. The characteristics of the jet 


*The value of the pressure equal to unity can be taken to be any pressure; experiments were carried out on mul- 
tiply-divide units built from the described schematic with P, = 1 atmos. 

**In the experimental model the unit containing the movable parts was mounted on a ribbon suspension to en- 
sure displacement of the nozzles without skewing; the force was transmitted, as is shown in Fig. 8, by membranes, 
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depend on the distribution of velocities in the output section of the nozzle. However, the effect of the above- 
mentioned factors is usually relatively small, and for free, turbulent jets,one can make use of certain average 
data on the basis of which one can calculate the characteristics of such jets. 


Additional processing of the experimental characteristics listed in earlier works dealing with the theory of 
free jets was carried out to clarify the numerical values of quantities determining the jet’s profile. The data ob- 
tained in this way was compared with the data of experiments carried out with nozzle-pipe elements in the course 
of the present investigation. The following considerations were kept in mind in the course of this. 


Fig. 9. 


1, In order to clarify the operating conditions of nozzle-pipe elements it is interesting to investigate the 
distribution of pressures in the jet. Works connected with the investigation of free jets, however, usually present 
characteristics of velocity distribution which, however, are also plotted on the basis of pressure measurements. 
Since the total pressure at any point in a jet, measured in gauge units, appears together with that velocity pres- 
sure, it is proportional to the square of the speed at a given point. Consequently, the practical boundaries of the 
jet, determined, for example, from the condition of equality on the boundary of the speed or pressure 0.5% from 
their maximum values, are different in those cases when velocity curves and pressure curves are being examined 
(compare, respectively, the quantities b and c in Fig. 9b; in Fig. 9b sw is the value of s for which W is equal to 
0.5% of Wa, in the given section of the jet). 


2. In conducting experiments with free jets the axis of the tube is usually placed parallel to the axis of the 
jet. Hence, although the velocity vectors in the jet are not collinear, in making use of formulas for the calcula- 
tion of the distribution of velocities, derived on the basis of experimental characteristics of pressure distribution, 
one can assume that at every point in the jet the axial component of velocity is equal to the total velocity. 


Let us first indicate the mean, calculated data with which the experimental characteristics of free jets and 
nozzle-pipe element: were compared. This data was obtained on the basis of fundamental relationships of the 
theory of turbulent, free jets, and the preliminary processing of experimental data on free jets.* The following 
numerical values were taken for the quantities determining the boundaries of the jet (see Fig. 9a): a= 8°40" 
(for the speed curve) or ¥_a = 7° (for the pressure curve); hzy = 5d which corresponds to 18 = 5°43’. In conjunc- 
tion with this it was assumed that the width of the boundary layer, determined as shown in Fig. 9a, is equal to 


b = hig + hig 8°40’ = 0,252h for 4 <5, (1) 


*Data on free jets obtained by Tryupel', Fertman, and from experiments at the Goettingen laboratory of aero- 
dynamics was processed further; for references to the respective sources see [4]. 
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b= 4 + hig840’ = $40,152h for (2) 


The following formulas were adopted for calculating pressures: 


Fig whereg > 5. (4) 


In these equations Pj, is the pressure before the nozzle, P ‘s the pressure at a given point in the jet, d is 
the diameter of the nozzle or its width (in the case of a slit-shaped nozzle), b is the width of the boundary layer, 
's is the coordinate of the point under consideration (Fig. 9a), h is the distance between the nozzle and pipe. 


Keeping in mind the characteristics of pressure and velocity distributions in sections of the jet normal to 
its axis, which are examined below we introduce the designations Wg, and Pg, corresponding, respectively, to 
the velocity and pressure at the center of the jet’s cross section, i.e., at the point in the section lying on the jet's 
axis, 


We can substitute values of b from (1) and (2), res vely, into (3) and (4). Thus, it is not hard to obtain 
the relationship P/P,x = f(s/b) and, since W/W, = the characteristic W/ Wax = f(s/d). 


Characteristics plotted from the equations presented above are shown in Figs. 10 and 11 as dotted lines. 
Figure 10a shows the jet’s profile, Fig. 10b shows the characteristic W/W, = f(s/b), Fig. 11 shows the character- 
istics P/Pg, = f(s/c) and P/Piy, = f(h/d). Experimental data for free, turbulent jets, and data obtained in the in- 
vestigation of nozzle-pipe elements are presented in these same figures. 


_ The following designations have been adopted in Fig. 10a: 1 points corresponding to Tryupel's data on 
a free jet emerging from a nozzle 90 mm in diameter (Re = 550,000, M = 0.26); 2— points corresponding to data 
from the Goettingen laboratory of aerodynamics, for a free jet emerging from a nozzle 137 mm in diameter (Re = 
= 380,000, M = 0.12); 3 — points corresponding to Fertman's data for a slit nozzle having a slit width of 30 mm 
and a slit length of 650 mm, corresponding to the case of the outflow of a free jet with Re = 70,000 and M = 
= 0.09; 4-—points obtained during the experimental investigation of a nozzle-pipe element with a nozzle hay- 
ing a circular cross section with diameter d = 0.8 mm, and a pipe with internal and external diameters, respec- 
tively, equal to 0.27 and 0.5 mm and length equal to 7 mm (Pip = 0.2-1 atmos+#M= 0.4-1, Re = 7,500-31,000); 
5 — points obtained by experimental means for a nozzle-pipe element having a slit nozzle with slit width 0.31mm 
and slit length 1.2 mm (M # 1, Re = 10,000); the pipe in this last case had the same dimensions as the one in the 
nozzle-pipe element with the circular nozzle described above. 


Together with the calculated characteristic W/W x = f(s/b) shown in Fig, 10b is the characteristic W/ Wax = 
= f(s/b), plotted on the basis of processed experimental data from the Goettingen laboratory of aerodynamics ob- 
tained for a nozzle 137 mm in diameter. The points designed in Fig. 10b by the numbers 1, 2 and 3 relate, re- 
spectively, to the values h = 0.6 m and h/d = 6.7, h = 0.8 m and h/d = 8.9, h = 1 m and h/d = 11.1.* Figure 11 
shows, besides the calculated characteristics, the characteristics 


P/Pax = f(s/c) and P/Pip = f(h/d), 


plotted from experimental data on nozzle-pipe elements the basic parameters of which were indicated above. 
The characteristics in Figs. 11a and b are plotted for a circular nozzle, the characteristics in Fig. 1lc and d were 
plotted for a nozzle having an exit aperture in the shape of a slit. In Fig. 11c are plotted points, for h/d = 6, ob- 
tained from the pipes crossing of the jet at three different transverse sections, shown in the nozzle's sketch. 


*Pressure in excess of atmospheric pressure. 
**For sections of a jet at a large distance from the nozzle, the characteristic W/ Wax = (f)(s/b), for low outflow 
velocities, approaches a straight line in the region between s/b = 0.3 and s/b = 1. 
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On the basis of the data presented in Figs. 10 and 11 we can conclude that the experimental characteristics 
of nozzle-pipe elements coincide sufficiently well with the experimental and calculated characteristics corres - 
ponding to the case of the flow of a free, turbulent jet. The design equations presented above are used subse- 
sequently in section 3, in the analysis of the variations in the pressure received by a pipe as it is moved in dif- 
ferent sections of the jet. However, it should be kept in mind that the conclusion arrived at concerning the co- 
incidence of the characteristics of nozzle-pipe elements and the characteristics of a free jet is valid only under 
the condition that the ratio of the area of the nozzle's output section to the area of the jet's section does not ex- 
ceed the values for which the nozzle-pipe elements described above were tested. The question of the effect of 


this ratio on the characteristics of pressure variations in nozzle-pipe elements is investigated separately in sec- 
tion 4, 


3. Characteristics of Pressure Variations Obtained during the Displacement of the 
Pipe Along Various Sections of the Jet "; 


Displacement of the pipe along the jet’s axis. When the apertures of the pipe and a nozzle having a round 
section are coaxial, the characteristic P = f(h) is defined for the main region of the jet (h/d < 5) by the relation- 
ship 


(0.3 + 0.44 (5) 


and for the initial region of the jet (h/d < 5) by the relationship 


P/Pjn = 1. (8) 


Relationships (5) and (6) are obtained from Eqs. (3) and (4) if we set s = 0 in the latter. The characteristic 
P/Pin = f(h/d), defined by expressions (5) and (6), is shown in Fig. 12a. From this characteristic, plotted in rela- 
tive coordinates, various families of curves P = f(h) can be obtained, For a given value of the nozzle diameter 
d, variations in Pj, result in changes in the ordinates of the characteristic P = f(h), proportional to the value of 
Pip (Fig. 12b).* However, variations in the nozzle diameter d for a constant value of P;,, is equivalent to varia- 
tions in the scale of quantities plotted along the h axis: for a given Pjp, the curves obtained for nozzles of dif- 
ferent diameters are shown in Fig. 12d. A row of other families of characteristics P = f(h) can be obtained if Py, 
and d are varied simultaneously, with their variations connected by some functional relationship. A feature of 
the characteristics shown in Fig. 12b and d are the relatively large values of the displacements h, corresponding 
to the working range of variations in P. These displacements have magnitudes of several millimeters or even 
tens of millimeters. What is significant is that the linear approximation of these characteristics is possible over 
sections several millimeters in length. 


Displacement of the pipe along the diameter of the section perpendicular to the axis of the jet. Substitut- 
ing the expression for b which follow from (1) and (2) into, respectively, (3) and (4) we obtain for the main re- 
gion of the jet (h/d > 5) the relationship 


s/d 
and for its initial region 
s/a\'h 


in 
For each given value of h/d the relationship P/Pjy, = f(s/d) is determined from (7) or (8). Figure 12e shows, 


*Taking the direction in which the tube approaches the nozzle to be the positive direction for changes in h, and 
choosing in a corresponding manner the origin of the h scale, we will obtain the characteristics shown in Fig. 12c. 
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Fig. 12. 


as examples, the characteristic P/Pj, = f(s/d) plotted for h/d = 1, 3, 6, and 10, Even by using any one of these 
characteristics a series of families of characteristics P = f(s) can be obtained. For example, for a nozzle with 
d = 0.8 mm and h/d = 6; i.e., with h = 4.8 mm, we will obtain, taking different values of Pip, a family of curves 

P = f(s) shown in Fig. 12f. Keeping Pj, = 1 atmos+ we obtain, for different d and for a given value of h/d, the e 
family of curves shown in Fig. 12g. For proportional changes in the quantities d and Pjp we obtain from the 
original characteristic P/Pj,, = f(s/d) with h/d = 6, shown in Fig, 12e, the family of curves P = f(s) shown in Fig. 
12h, etc. 


Movement of the pipe in a plane forming an angle y with the jet’s axis. We will limit our investigation 
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to the case where the pipe moves along the major axis of the ellipse obtained when the main region of the jet is 
cut by a plane forming an angle y with the jet’s axis (Fig. 121), Moreover 


(s,/ d) sin @ 
1— 2.85 (; + 0.304 | (9) 


0.3 + 0.14 (hy (8, / d) cos 9) 


The values of the quantities y, sy, hy are given in Fig. 121; the rest of the quantities in (9) are the same 
as those appearing in earlier expressions. Figure 12i also shows the limiting values of sy, being, respectively, 
Sy ,max 4d 8% max. When 0 < 8» < sy,max the sign in front of the term (sg/d)cos¢ in (9) should be taken as 
plus, while when 0 < sy < Sy,max, the sign should be taken as minus. The quantities Sy,max 4Nd 8, max can 
be determined from a drawing of the jet's profile (Fig. 10a) or they can be calculated from the equations 
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cos @ (tg @ + 0.152) * 
The last expression in (10) is valid only in those cases when the plane of the jet's section does not touch 


and may 


(10) 


the jet's initial region, which corresponds to the condition tang = wey When tang < By one 
should use the equation 
 @ (11) 


The derivation of expressions (9), (10) and (11) is given in Part 1 of the Appendix. We note that when ¢ = 
= 1/2, we obtain Eq. (7) as a particular case of (9); here hy = h and sy = s, When g = 0, we obtain formula (5) 
as another particular case of (9); here hy 4 Sy = h. 


Figure 121 shows the characteristics P/Pip = f(sy/d), obtained from Eq. (9) for hg/d = 6 when ¢ = 0, 5, 10, 
30, and 90°, 


Further expansion of the class of characteristics P = f(s) obtained in sections of an axially symmetrical, 
turbulent jet. Families of characteristics P = f(s) different from the ones considered above can be obtained if 
the pipe is displaced along the chords rather than the diameter of a section perpendicular to the jet’s axis, In . 
this case, the characteristic Pjy/P = f(s/d) is determined, for a given h/d, from the formula 


2 
_ 2.83 + Dina 
1+ 0,304h/d 

0414h/a 


(12) 


Pin 


in which s, is the distance from the section's axis, measured along the chord, and a is the distance of the given 
chord from the jet's axis. For sy max/d we obtain, in this case, the equation 


The pipe can be displaced along an axis parallel to the jet's axis and separated from it by a distance s. 
The family of characteristics obtained in this case is described by Eq. (7); one must only take in this equation 
the value of s/d as constant and consider h/d as the variable quantity. The minimum value of the ratio h/d for 
which the pipe, moving in the direction of the nozzle, reaches the boundary of the jet, is equal to 


s 
0.5 
h d (14) 
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The derivation of Eqs, (12)-(14), as well as further comments on the use of Eq. (7) in the case under con- 
sideration, are given in Part 2 of the Appendix. The class of families of characteristics P = f(s) can be widened 
even further if the pipe is displaced along the chords of curves formed by the intersection of inclined planes with 
the jet. 


All the relationships presented above refer to the case where the pipe's axis is parallel to the jet's axis 
(Fig. 12j) where, as it was mentioned, it is assumed that the total velocity pressure of the jet acts on the pipe. 
The pipe can be situated at an angle @ to the jet's axis (Fig. 12k). However, for large angles of inclination of 
the pipe, the disturbances introduced into the air flow by the pipe become appreciable, and in performing cal- 
culations it is necessary to introduce an additional correction to the value of the velocity head, determined from 
the speed Wg = Weosg- In such a case, for example, (7) is transformed into the expression 


r s/d 
| 1 — 2.83 a) | 


Pin =& 0.34 0.14h/d | (18) 


where k is a correction coefficient. On the basis of experiments conducted with a pipe having internal and ex- 
ternal diameters equal to, respectively, 0.27 and 0.5 mm (nozzle diameter 0.8 mm) this coefficient can be de- 
termined from the expression* 


k = (16) 


In this way, another variable quantity 6 is introduced into the expression for the functional dependence of 
P/Pj,, on s/d and h/d determined previously from Eq. (7). In a particular case, when @ = 90 — g, in (15) cos @ is 
replaced by sing. 


Expressions (5), (7), (9), and all subsequent ones were obtained for the main region of an axially symmet- 
rical jet. Analogous expressions can be derived also for the initial region of the jet. 


Calculations of the parameters of a nozzle-pipe element. We will limit ourselves to examining the char- 
acteristics obtained when the receiving aperture of the pipe is moved along the jet's axis, and determined by 
Eq. (5).** First let us present formulas for the calculation of the distance from the nozzle to the pipe hy and the 
pressure before the nozzle Pj, for which the curve P = f(h) passes through the given points A (P9,hg) and B(Py,hy + 
+ Ahy) (Fig. 12m)***; 


17 
V 


(18) 


Let us illustrate the use of Eqs. (17) and (18) by the following numerical examples. Let Py = 0.65 atmos+, 
P,=0.38 atmos+, Ahy = 2 mm. Let us take d= 0.8 mm. From Eqs. (17) and (18) we find hy » 4.7 mm and Pjp = 
= 0.83 atmos+, 


We can specify that the characteristic P = f(h) passes through the points A (Po,hg), B(Py,hg + Ahy), C(P2,hg + 
+ Ahy), ..- , M(Pyp,hg + Ahm), N(Pn,ho + Ahn), the coordinates of which satisfy the condition (19). 


We note further that if it is more convenient, from a practical point of view, to specify the characteristic 
P = f(h) by the angular coefficient (dP/dh), of the tangent to the characteristic at the point P = Pp, then the quan- 
tities hy and Pj, can be calculated from Eqs. (20) and (21). 


*For more detail see next section. 

**Only in this case does the calculation of the nozzle-pipe element's parameters become sufficiently simple. 
In other cases the calculation reduces to the solution of nonlinear equations, which can be carried out only by 
graphical means. 

***The derivation of formulas (17), (18), and formulas (19)-(21) which are given further, is given in Part 3 of 
the Appendix. 
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Comments about the characteristics of slit-shaped nozzles. In part 2 it was shown that the characteristics 
of slit nozzles (Fig. 3a) obtained when the receiving aperture of the pipe is moved along the jet’s axis, or those 
obtained when the pipe moves in a transverse section of the jet along the straight line parallel to the z axis 
(Fig. 11c) are also closely approximated by the design equations presented earlier. In these cases d, in the given 
equations, is the width of the slit. For practical applications it may be important to determine the characteris- 
tics obtained when the pipe's receiving end is moved along the major axis of the jet’s transverse cross section 
(along the axis of the jet's transverse cross section, lying in the xy plane, see Fig. lic). By the use of just such 
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Fig. 13. 


characteristics, the universal nonlinear converter, using nozzle-pipe elements, described in part 1 was designed. 
Experimental investigations showed that the distribution of pressures in the boundary layer of the jet flowing out 
from the slit nozzle (we have in mind the inclined portions of the characteristics), differs little for the cases 
where the pipe moves along the transverse cross section's central axes parallel to the z axis and the x axis. 
Figure 13a shows, as an example, the corresponding regions of the characteristics P = f(s,) and P = f(x) plotted 
for an element with a nozzle having a cross section 0,31 mm wide and 1,2 mm long, when h = 1.8 mm and Py, = 
=1 atmos+.Figure 13b shows the experimental characteristics P = f(s,) for the same element with h = 0.3 mm and 
Pin = 1, 0.6, and 0.2 atmos+. These characteristics are very close in their shape to the theoretical characteristics 
shown earlier in Fig. 3. The experiments were carried out with a pipe having an internal cross section diameter 
of 0.27 mm and an external diameter cf 0.5 mm. 


4. Investigation of the Deviations of Nozzle-Pipe Element Characteristics from 
the Characteristics of Free, Turbulent Jets 


Deviations of the characteristics of nozzle-pipe elements from the characteristics of a free jet can be 
brought about by the following factors, With relatively large dimensions of the pipe's receiving aperture (with 
small nozzle dimensions) the presence of gradients in the velocity head pressures in the jet’s section assume signi - 
ficance, In this case the pressure measured by the pipe is some mean of the pressures acting on the plane of its 


receiving aperture. Moreover, when a tube of relatively large dimensions {is introduced into the jet, the variation 
in the pattern of flow in the jet in front of the pipe, which is brought about by the action of aerodynamic factors, 
can have an effect on the pressure measured by the pipe. 


Investigation of the effect of the pipe’s receiving aperture. First we will assume that in the region in front 
of the pipe’s receiving aperture the velocity field is the same as it would be in a free jet, i.e., that the pipe does 
not produce distucbances in the jet, and we will assume that the process by which the pipe slows down the air is 
isentropic. (Lier, possible deviations from this idealized case will be discussed.) Under the above assumptions 
the energy of the slowed-down flow is characterized by the quantity = | Pdf, B P,, f+ where Pye is the pressure 


2 
Fig. 14. 


For positions of the pipe in sections of the main region of an axially symmetrical jet shown, respectively, 
in Fig. 14a, b, c, and d, the quantity Pye can be calculated from the equations given below, the derivation of 
which is gives in Part 4 of the Appendix. The following notation has been introduced in Fig. 14: b — radius of 
the jet*s section, a ~ distance from pipe’s axis to jet's axis, s and g — polar coordinates of points in the jet, which 
lie in the field of the pipe’s receiving aperture. Also, in the formulas given below there appear, as before, d — 
nozzle diameter, h — distance between the ends of the nozzle and pipe, Pin — air pressure supplied to the nozzle. 


For the case given in Fig. 14a: 


2 s/d 


Equation (22) also holds for the case shown in Fig. 14b with the sole difference that the upper and lower 
limits of integration are replaced, respectively, by b/d and a/d—r/d. For the location of the pipe in the jet, 
shown in Fig. 14c: 


r—a 
2d?P, 2 sid 
Pye = T+ 0.304 hid Ty's 
‘ 0.340.14 h/d 
0 


*We will note once more that Pye and P are the pressures of the velocity head and are, at the same time, the 
total pressures measured in excess units, 
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Equation (23) can also be used for the case shown in Fig. 14d if the upper limit in the first integral on its 
right side is taken as r/d and the second integral is set to zero. 


For the location of the pipe in the jet shown in Fig. 14a, calculations from Eq. (22) are carried out for 
given values of h/d and r/d in the following way. The limits of variations in a/d are taken into account (in the 
case being considered r/d < a/d< (b—1)/d and a few values of a/d which lie within these limits are selected, 
The limits of integration are calculated for every chosen value of a/d, and several values of s/d which fit be- 
tween these limits are taken. The subintegral function is calculated for them. Its graph is plotted point-by~point, 
and a planimeter is used to determine the value of the integral in Eq. (22). These calculations and plots are re- 
peated for other chosen values of a/d. The considerations presented above also apply to the calculation of the 
second integral on the right side of (23). The subintegral function of the first integral on the right side of (23) 
does not contain a/d and the calculation of this integral reduces to merely picking several values of s/d between 

Plo. the limits 0 < s/d< (r/d —a/d), plotting the curve of the sub- 
w integral function, and determining the area bounded by the curve 
by a planimeter, 


\ The method described above was used to calculate the 
as / characteristics P = f(h) obtained for displacements of the pipe 


along the jet's axis, for nozzle-pipe elements with nozzle di- 
ameter d = 0.8 mm and diameter of the pipe’s internal section 

ata i 50 dp = 0.27 and 0.8 mm. These characteristics were compared 
y | with the characteristic P = f(h), calculated from Eq. (5) under 


the assumption that dp = 0. All three characteristics are shown 
in Fig. 14f. 


Given experimental methods for determining the effect of 
Fig. 15. the pipe dimension on the characteristics of a nozzle-pipe ele- 
ment. Considerations regarding the effect of the disturbances 
introduced by the pipe into the jet flow. For a nozzle-pipe ele- 
ment with a pipe having dp = 0.27 mm, Fig. 14g shows the experimentally determined (solid line) characteristic 
P/Pin = f(h/d) and the characteristic P/P;,, = f(h/d) calculated by the method described above (dotted curve). 
These characteristics are practically the same. Figure 14h shows the experimental (solid) and calculated (dotted) 
characteristics P/Pj, = f(h/d) for a nozzle-pipe element with dp = 0.8 mm (as before, the nozzle diameter d = 
= 0.8 mm). This experimental characteristic passes lower than the experimental curve for the pipe with dp = 
= 0.27 mm (Fig. 14g). Thus, the tendency for further decrease in the value of P with increasing dp is confirmed, 
but the observed values of P are higher than the calculated ones. The latter can be explained by the fact that 
beside the effect of the area of the pipe's aperture cross section on the flow characteristics of nozzle-pipe ele- 
ments, which was discussed above, the effect of disturbances introduced by the pipe into the stream can, as has 
already been mentioned, begin to appear. When the losses during deceleration are taken into account, the pressure 
received by the pipe must be smaller than that obtained from calculations carried out for the case of isentropic 
deceleration, However, an increase in pressure can occur as a result of the deformation of the streamlines in the 
jet (as a result of which the energy of a greater part of the jet than that which was taken into account in the de- 
scription of the calculating system, is transformed into the velocity head pressure. An increase in pressure can 
also occur as a result of the introduction into the deceleration region under consideration, of particles from layers 
adjacent to it (Fig. 14e). In this case, it is probable that not only the internal but also the external pipe diameter 
affects the characteristics of the nozzle-pipe elements.* For the pipe with dp = 0.27 mm, the outside diameter 


*The pipe’s length and the profile of its leading end can also affect its streamline conditions, All experiments 
were conducted with pipes of length 2 = 7 mm, and having their end surfaces perpendicular to the pipe axis. 


was equal to 0.5 mm for the experiments described above, while for the pipe with dp = 0.8 mm the outside di- 
ameter was 1.26 mm. The area corresponding to the outside diameter for the first of these pipes was 2.5 times 
smaller, and for the second 2.5 times larger than the cross-sectional area of the nozzle. Experiments were car- 
ried out with a pipe having an external diameter equal to 3 mm and dp = 0.8 mm to determine what effect fur- 
ther increases in the pipe’s external diameter have. The obtained characteristics (1 — for a pipe with a 3 mm 
external diameter, 2 — for a pipe with a 1.26 mm external diameter) are shown in Fig. 14i; they practically do 
not differ from each other. 


All the experimental characteristics described earlier were obtained with the receiving aperture of the pipe 
situated in a plane perpendicular to the jet's axis. As it was already noted, if the pipe’s angular position is varied 
it is necessary to take into account in the calculations the variation of the pipe’s streamline conditions as a func - 
tion of the angle © (Fig. 12k). Figure 15 shows the experimental characteristic P/Pg -» = f(@), obtained by rotat- 
ing the pipe about a vertical axis passing through the center of its receiving aperture, with this point situated on 
the jet*s axis. Experiments were carried out with a nozzle having d = 0.8 mm, and a pipe with internal diameter 
0.27 mm and external diameter 0.5 mm, with h/d = 6. The dotted line in Fig. 15 shows the corresponding cal- 
culated curve obtained by calculating k from Eq. (16). 


APPENDIX 


Derivation of Formulas Used in the Calculation of the Characteristics of Nozzle- 
Pipe Elements 


1. Derivation of expressions (9)-(11). Let h = hg + Ah (Fig. 16a). Then from AABC: s = sysing, Ah = 
= s_ cosy and h = hy + sycosy. .rom ADEF: b ~—d/2 = (hg + Ah) tan(a/2) and, hence, b = (hg + Ah) tan(a/2) + 
+ d/2 = (hg + Sg cosy) 0,152 + d/2. If h = hy — Ah (Fig. 16b) then from AAB'C’ it follows, as before, that: s = 
= sing and Ah = sg cosy. However, in this case h = hy —s,cosy. From AD*E'F* we find b = (hy —syeosg) x 
x (0,152) + Yod. Making use of Eq. (7) for determining the pressure, respectively, at points C and C’, and 
substituting into (7) the values of s, h, and b shown above we end up with (9). 


Expressions (10) are obtained thus, From AAGH (Fig. 16c): by = Ahymax tang. From AGEF: c = (hy + 
+ Ahymax) tan(a/2) = (ig + Ahymax) 9.152. In this case by = 0.152 (hg + Ahymax) + 4/2. Eliminating Ahymax 
from the two expressions for by we find by = (0.152h, + d/2) tan g/(tan g — 0.152). From AAGH: S9smax = b,/sin ¢. 


Substituting the value of b, into the last expression we obtain the first equation in (10). An equation for 
determining Sy, .,,,,/d is derived in an analogous manner. However, the last is valid only up to the limiting 
position of the section plane shown in Fig. 16d, where this plane intersects the jet's boundary at the point with 
abscissa 5d, which corresponds to the transition from the initial region to the main region of the jet. For this 
position we find from AOKL (Fig. 16d) f = 5dtana/2 = 0.76d. In this case NL = 1.26d and from AALN it follows 
that tang = 1.26/(hg/d —5). When tang < 1.26/(hg/d—5), which follows from Fig. 16e, 0, = (Ne — 5d)/cosg 
and $95 max/4 = (g/d — 5)/cos ¢, i.e., we obtain expression (11). 


2. Directions for the derivation of Eqs. (12)-(14) and the use of Eq. (7) in the case of the displacement of 
the plpe along the axis parallel to the jet's axis. Equation (12) is obtained from (7) by replacing s by +a 
(Fig. 16f). Expression (13) follows from the relationship sx max = Vb" —a* (Fig. 16g). An examination of Fig. 16h 
leads us co conclude that it is possible to make use of Eq. (7) in the case of the displacement of the pipe along 
an axis displaced by a distance s from the jet's axis. From AEFD we find (s — ¥gd)/hyyj,, = tana /2 = 0,152, 
whence (14) follows. 


3. — of Eqs. (17)-(21). First we will indicate the method of eisinis Eqs. (17) and (18). From (5) 
1 0.3 + 0.14 [hy + Ahyj/d 
Po = Pin and Py = Pin . Hence, Po/P; + 
(0.3 + (0.3 Ah 0.3 + 0.14h,/d 
Solving the last equation with respect to hy,we arrive at (17). Substituting the obtained value of hy into the ex- 


1 
1 
pression Pg = Pin (0.3 + 0.141 We (1S). Let us further examine the characteristic 2n. 


Analogously to the way in which Eq. (17) was derived we can find the expressio Aha /( ¥Po/Ps foes 


Equations for each section up to hy + Ah, can be treated analogously. As a result we will obtain conditions (19). 
Equation (20) is derived as follows: we differentiate (5) from the expressions 
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0,28 P 4 
we eliminate Pin and, taking into account that (dP/dh)y < 0 and, therefore, | (dP/dh)»| = —(dP/dh)p, we arrive at 
(20). Equation (21) is obtained by substituting the vale of hy from (20) into the expression Py = Pin 
4. Derivation of Eqs. (22) and (23). When the tube is situated in a section of the jet, as shown in Fig. 14a, 


the pressure is the same at all points in the shaded area [it is determined from Eq. ('7)], but the area of this re- 
gion is equal to f = 2sy ds, Because cos ¢ = (a* + s* —r*)/2as, then f = 2s cos([a® + s* —1r*]/2as) ds. 


Fig. 16. 


An elementary force acting on this area is equal to 


P 034 0,i4h/a 


Replacing s, a and r by the ratios s/d, a/d and r/d, integrating with respect to s/d between the limits of 
a/d—r/d and a/d + r/d,and dividing the resultant force thus obtained by the area of the pipe’s aperture rr*, we 
obtain (22). It follows from Fig. 14b that the derivation presented above includes the gtven case, s/d = a/d —1/d 
and s/d = b/d are only the limits of integration. 


For the case given in Fig. 14c, one takes the sum of two integrals, one of which corresponds to the interval 
0<s<-14r< a, and the other to the interval r —a < s < r + a (the shaded region in Fig. 14c corresponds to the last 
interval), The derivation of Eq. (23) is analogous to the derivation of Eq. (22). The second integral in Eq. (23) 
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is the same as the integral in Eq. (22), and it corresponds to the second of the intervals indicated above. The 
elementary areas for the first interval are determined by the expression f = 2rsds, which is used in determining 
the first integral in the right side of (23) instead of the expression f = Qsarc cos ([a® + s* — r*}/2as) ds, which was 
given before. The coaxial location of the apertures of the pipe and nozzle (Fig. 14d) can be considered to be a 
particular case of the relative location of the pipe in a section of the jet shown in Fig. 14c. In this case, only 
the first of the integrals in Eq. (23) remains and the upper limit of integration is equal to r/d, since a = 0. 
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DESIGN PRINCIPLES FOR CONTACTLESS REMOTE-CONTROL SYSTEMS 
WITH EXPONENTIAL CONVERTERS*# 


V.A. 


(Moscow) 


Design principles for proposed, new, signal time-separation systems of remote 
control and remote signalling are described. Formulas for determining the main 
parameters of such systems are presented. 


The prospects of replacing remote-control systems with electromagnetic relays by the more reliable con- 
tactless systems withtime- or frequency-separation of signals have become very good recently. 


One of the basic elements of contactless, signal time-separation systems of remote control consists of a 
multiwinding transformer with a ferromagnetic core having a rectangular hysteresis loop, crystal diodes, and 


capacitors [1]. With relatively slow operation, difficulties often arise in such a system in matching it to a com- 
munication channel. 


In signal frequency-separation systems of remote control the simplest and most promising frequency dis- 
criminator is the electric circuit (or band filter) containing an inductance and a capacitor [2]. 


From the point of view of complexity the indicated contactless elements 


of systems withtime- and frequency-separation, for the same requirements, are 
approximately equivalent. 


This paper presents the design principles for signal time-separation systems 
of remote control and remote signalling, the basic element of which is a new, 
simple, exponential converter which has already found application in telemetry 
systems [3]. The exponential converter, consisting of resistors, a capacitor, and 
a diode, permits one to design, in a number of cases, simpler systems with time- 
and code-separation of signals, without the use of inductances or transformers. 
Such systems can be well matched to a communication channel, they are very 
stable and can be made smaller and lighter. This makes their prospects of de- 


velopment very promising. 


The simplest circuit of an exponential converter which transforms the magnitude of a resistance or capaci- 
tance into a pulse duration is shown in Fig. 1. The converter consists of a bridge, one arm of which contains re- 
sistors Ry and Ry, while the other contains a resistance R and a capacitor C. A diode is connected across a diagonal 
of the bridge. The forward resistance of the diode Rp and the resistance R are chosen from the conditions Rr « R 
and Rg « R. The diode's reverse resistance R,¢ must be considerably higher than the resistance R. 


When rectangular pulses (Fig. 2a) derived from a blocking-oscillator or some other pulse generator, are 
applied to the bridge the voltage across capacitor C changes as shown in Fig. 2b while the current through the 
capacitor changes as shown in Fig. 2c. At time t, the voltage across the capacitor becomes equal to the voltage 


Fig. 1. 


*The author presented a précis of this article at the All-Moscow Seminar on Automation and Remote Control 
held at the Inst. Autom, and Remote Control, June 11, 1958. 
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Fig. 2. 


across resistor Ra, the diode begins to conduct and, for the chosen 
parameters and Ry « R, the capacitor stops charging. The duration 
of the converter’s output pulses is determined by the expression 


te= RCN (1) 


The trailing edge of the pulses lengthens if the diode does not 
have an ideal piecewise linear characteristic. In practical converters 
durations of the trailing edge of the order of 1% of the maximum dura - 
tion of the output pulses can be obtained. The duration of the trailing 
edge will be 


at = Fac (2) 


where E is the voltage amplitude of the primary pulses, AE are ap- 
proximate voltage variations along the diode’s characteristic, for 
which the diode's internal resistance changes from R/2 to 2R. 


In the simplest system of remote control which makes use of 
a distributing code with n commands, each command is given a dif- 
ferent pulse duration by the exponential converter. The object switch 


K in the transmitter (Fig. 3, left) switches resistors and gives the pulses a particular duration t,j (resistor R or 
capacitor C can be varied). The converter's output voltage is delivered from the load resistor R; « R to the 
linear repeating relay LR, which can be a transistorized, electronic relay, etc. 


The right side of Fig. 3 shows the circuit of a three-command receiver. The receiver consists of a dis- 
tributor in the form of a complex, exponential converter with four identical RC networks, a differentiating stage 
DS and a separate coincidence circuit CC for each command. 


The time plots of voltages across load resistors Ryy, Ri, Rug and Ry, are shown in Fig. 4b, c, d and e, while 
the output pulses from the distributor, which are delivered to the inputs of the coincidence circuits, are shown in 


Fig. 4f, g and h. 
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Fig. 3. 


A pulse corresponding to the trailing edge of the pulse transmitted from the transmitter is delivered from 
the output of the differentiating stage to all the second inputs of the coincidence circuits. If the duration of the 
pulse sent from the transmitter is greater than the duration of the pulses across the load resistor Ry, and smaller 
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than those across Ry», then the pulse from the differentiating circuit coincides with the pulse of the distributor's 
first channel, which corresponds to the transmission of the first command, etc. The coincidence stages can be 


circuits designed with diodes, transistors, etc. An output seal-in relay is connected, if required, to the output 
of each stage. 


The input pulses can be produced not only by a pulse generator in the transmitter (Fig. 3) but also by a 
pulse generator in the receiver, similarly to the way it is done in the VST~-1 telemetry system [3]. In many cases, 


for example, in the use of frequency channels, it is necessary to have generators of primary pulses, synchronized 
by well-known means, both in the transmitter and receiver. 
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Fig. 4. Fig. 5. 


The maximum number of transmitted commands with a distributing code is 


_ max 
n=K Btay —1, (3) 


where Atay = (Atmin + Atmax)/2, K is a factor of assurance (K » 0,7-0.8). The number of commands n, cal- 
culated from Eq. (3), can be increased several times if special forming circuits are introduced. 


It should be noted that multichannel, signal time-separation systems formed by the selection of different 
levels of trigger voltages, were repeatedly described in the literature, while the mechanical analog of such sys- 
tems was already realized in Prof. P.A. Molchanov's radiosonde. Thus, for example, in the circuit proposed by 
F.E. Temnikov, the separation of signals was achieved by the use of a special sawtooth generator, and a “bunch” 
of null-elements and rheostats connected in series [4]. 


However, such systems had relatively complex null-elements, auxiliary sawtooth generators, and other 
equipment with insufficiertly stable parameters. In contrast to this, circuits with exponential converters do not 
contain unstable components which significantly affect the system's operation. Thus, the use of exponential con- 
verters permits one to realize simple and very stable, remote-control systems, 


Figure 5 shows the circuit of a transmitter for a remote-control system with a decimal code, designed for 
the transmission of 100 commands, It contains two RC networks, two diodes and, correspondingly, two switches 
S; and S,, designed for setting up tens and units. For the transmission of, say, order No. 67 switch S, sets at posi~ 
tion 6 and switch S, sets at position 7. 


The output of the exponential converter is connected to a differentiating network CgRq for sending pulses 
into the line when the first and second diodes open. The principle of pulse separation here is the same as in a 
multichannel, telemetry system using exponential converters [3]. 
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The receiving unit of a remote-control system must have a distributor with individual outputs and coinct- 
dence circuits for each ten and unit, designed like, for example, the receiver in Fig. 3. 


The circuits of the transmitters shown in Figs. 3 and 5 can transmit during one cycle only one command. 
For transmitting two or more commands during one cycle a correspondingly larger number of RC networks is in- 
troduced in the transmitter, for example, as is shown in Fig. 5. 


SUMMARY 


1, Exponential converters permit the production of simple and very reliable, contactless systems of remote 
control and remote signalling with signal time division, without the use of components with ferromagnetic cores. 


2, Compared to other systems with time- and code-separation of signals, the circuit of the transmitter with 
exponential converters can be simplified considerably if the number of commands transmitted during one cycle 
is decreased. 


3. For the reception of one command in the systems being considered the receiver's distributor can have a 
simplified exponential converter with two diodes and one coincidence circuit. This simplifies considerably the 
receiver, as compared to other systems with signal time-separation, and opens up wide areas of application for 
exponential converters in remote-control systems for distributed objects. Frequency systems of remote control, 
which are widely used for the remote control of distributed objects [2] possess this same property. 


4, Pulse duration and period of cycle in an exponential converter can be chosen in the range of tens of 
microseconds to tens of seconds. This permits the variation of the system's response over a wide range, and sim- 
plifies matching the system to practical transmission chamnels. 


LITERATURE CITED 


{1] V.A. Zhozhikashvili and K.G. Mityushkin, “Concerning the operation of computing~switching circuits 
using magnetic elements with PPG,” Automation and Remote Control (USSR) 16, 4 (1955). 


(2] V.A. Din and K.P. Kurdyukov, “Concerning frequency methods for the remote control of distributed 
objects,” Automation and Remote Control (USSR) 19, 2 (1958).* 


(3] V.A. Din and A.L Novikov, “New design principles for telemetry systems with pulse-time and pulse- 
width modulation,” Automation and Remote Control (USSR) 19, 8 (1958).* 


[4] G.M, Zhdanoy, in the book: Telemetry, Part 2 [in Russian] (Gosénergoizdat, 1953). 


Received October 16, 1958 


*See C.B. translation. 


452 


t 
] 
t 


ANALYSIS OF THE OPERATION OF PULSED MAGNETIC ELEMENTS 
IN CONTACTLESS REMOTE-CONTROL DEVICES* 


I.V. Prangishvill 
(Moscow) 


This article provides an analysis of the operation of pulsed magnetic elements 
with rectangular hysteresis loops when the elements are fed by a sinusoidal voltage of 
industrial frequency. A series of characteristic relations is deduced. A method for 
calculating a single element and a pulse relay is given. 


Among the contactless devices of automation and remote control, more and more use is being made of 
pulsed elements, with magnetic cores that have rectangular hysteresis loops (RHL). 


Magnetic elements with RHL together with semiconductors serve as the basic elements of the contemporary 
contactless devices used in automation and remote control. They are used in the various loops of remote-control 
devices as pulse amplifiers or repeaters, memory elements, and so on. 


Magnetic pulse elements with a large amount of positive feedback possess the characteristics of contact~- 
less polarized relays [1]. Such elements are used as output and intermediate relays. 


It is possible to study two groups of circuits composed of pulsed magnetic elements. The first group con- 
sists of circuits which have series-connected loads, these are known as choke circuits [2]. The second group con~ 
sists of circuits which have parallel-connected loads, these are known as transformer-type circuits [3-6]. 


The elements of the second group have several advantages over those of the first group. First of all, there 
is the ability to witlstand large oscillations in the input voltage. It is because of this that the elements of the 
second group are used, in preference to those of the first group, in applications which take place during the de- 
velopnient of remote-control devices. It is for this reason, that the present article will be devoted only to the 
study of pulse circuits with transformer-connected loads, 


Of all the possible methods of feeding such elements, the simplest and most reliable seems to be that of 
using an alternating voltage of industrial frequency. In spite of this, this method of supply leads to a relatively 
large size of circuit element. The core size of the pulsed magnetic elements now in use depends mainly upon 
the frequency and form of the input pulses. 


When building large-capacity devices, where a decrease in the size of the elements to be used in the de- 
vice is significant, it is expedient to raise the frequency of the input pulses, and also to make use of special 
pulse shaping devices. However, for devices of relatively small capacity, for which the question of size is not 
of basic importance, the first consideration becomes that of achieving the maximum simplicity and reliability 
in the pulse supply. 


These considerations determine the choice of the type of impulse generator to be used with a device of 
relatively small capacity (20-30 objects) operating on industrial frequency ac. The half-waves of the sinusoidal 


*Presented at the General Moscow Seminar on Complex Remote-Control Systems held in 1956. 
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Fig. 1. Simplified circuit 
diagram of a magnetic ele- 
ment with sinusoidal input. 


ac voltage may be used as the input pulses for the device. The operation of 
elements possessing RHL under these conditions is characterized by definite 
specifications and it differs from the operating conditions that exist when the 
circuit is supplied by current pulses, as is the case in computer techniques.. 


When elements with RHL are fed from a half-wave alternating voltage, 
the stable operation of these elements requires that there be a sufficiently 
rapid and complete magnetic polarity reversal of the core magnetization in 
order that the development of the necessary signal power at the output ele- 
ments be insured. In connection with this, the time required for complete 
polarity reversal of the core magnetization ty), , may take on various values 
within the half-period of the supply voltage. The amplitudes and the effec - 
tive values of the voltage and power of the output signals depend upon the 


magnetization time. In its turn the time required for a complete magnetic polarity reversal is determined by 


the amplitude of the input voltage and by the parameters of the magnetic element. With a given magnitude of 
output signal, the length of the time interval that is required for complete magnetic polarity reversal of the core 
can be determined, and from this we can calculate all the parameters of the magnetic pulse element. There- 
fore, in order to analyze and calculate the values of these elements, it is necessary to determine the dependence 
of the output voltage and power upon the time required for a complete reversal of the magnetic polarity of the 
core when the elements are fed by an alternating sinusoidal voltage. 


g+ 


Fig. 2. Graph of the dependence of Uy max/U'y max: UH/U'H max: 
upon the magnetization time 


We will now turn our attention to the determination of this dependence. While doing this, we will approx- 
imate the actual hysteresis loop by a rectangular loop, and we will make use of the fact that the integral of the 
electromotive force existing in the winding taken over the time interval tmag required for a complete reversal 


454 


of the magnetic polarity of the core is a constant, and is independent of the magnetization time: 
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edt =—2B,WS, -10*, 


where e represents the electromotive force in the winding being examined, B, represents the residual induction, 
taken as being equal to the saturation value, W is the number of turns in the winding and S, is the transverse 
cross section of the core. 


This expression permits us to determine a series of relations characterizing the operation of magnetic ele- 
ments for various magnetization times, and allows us to simplify the procedures required to calculate these ele- 
ments, In order to determine the indicated relations let us examine the simplified circuit of a magnetic element 
with RHL, as shown in Fig, 1. We will neglect the resistance of the winding and the rectifiers. 


L The determination of Upj max = f1(tmag), Where Uy max tepresents the maximum value of the voltage 
(or the electromotive force) in the load windings when the core is not saturated. 


Applying Eq. (1) to the load winding, we ge* 


dt = —2B,W,S, -107*. (1a) 
0 


On the other hand, since the input voltage is sinusoidal we get 
ag 


Uy, dt = — Tamar (4 — €08 Wtmag)- 


If we let 2B,Wy;S 4° 10~* = m, we may write 


Un.max = (2) 


Let us designate by U"}zmax the value of Uj max for the case where trag = T/2. Then, from (2), we get 


max = 9-5 mw. (2a) 
Dividing (2) by (2a) we find: 
Uy max “i 2 : (3) 
i -— cos oot na g 


IL The determination of Uy; = fg(tmag), Where Uy is the effective value of the electromotive force upon 
the load: 


U al 4 dt =U max Bin 
Substituting the value of Uy; 4x obtained in (2), we obtain 


ML The determination of uyy = f3(tmag)» where uy; is the instantaneous value of the load voltage at the 
end of the magnetization time, i.e., for t = thag. 


mo i 
Un = 7 8in 2otmap 
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Using the fact that uy; = Upymax in wtmag, in Eq. (3) we obtain the expression 


Uy 2sin 


; = 5 


IV. The determination of Pj; = f4(tmag), Where Py is the average load power. 
Using Eq. (4) we obtain 


(me)? 


Pu= R uf (1 — cos of, 


(Spas sin 20tmag (8) 


Let us designate by P’}; the value of Py; when tag = T/2. Then from (6) we obtain 


eR, (6a) 


P, = 


Dividing (6) by (6a) we get 


t 


(6b) 


Graphs of the expressions (3), (5) and (6b) are shown in Fig. 2. 


With the aid of expressions (1)-(6) it is easy to analyse and to simplify the calculation of magnetic pulse 
elements under various operating conditions. 


The basic element of contactless remote-control devices consists of the circuit of a single element with 
resistive load and without feedback paths (Fig. 3). Therefore, before going to the calculation of pulse relays, 
we must first stop and calculate the operation of this type of element. While we are doing this, let us agree to 
call the magnetization of the core due to the input impulse the “operation” of the element, and the time re- 
quired for this we shall designate as tprep. 


Methods of Calculation of a Magnetic Pulsed Element with Resistive Load 


We shall consider the reverse resistance of all rectifiers as being infinite. We shall only take account of 
the forward resistance of a rectifier when it occurs in the preparatory circuit, where its magnitude has a signifi- 
cant value. 


For the load circuit, from Fig. 3, 


Uy = Eq (7) 
where E;;, Uz; and ky represent the effective values of the electromotive force, voltage and current. 
The voltage drop in the load windings may be represented by 


l 
In Ry Age = 170-10 Aly Wa, (a) 


where A is the current density in the load winding in amp/mm’*, 1, is the average length of one turn of the load 
winding in cm and 170- 10~° is the specific resistance of copper at 20°C in ohm-mm*/cm. 


Substituting for the value of kyRy}; in (7) the expression found in (7a), and for the value of Ey; in (7) the 
expression found in (4), we obtain 


Solving the last equation for W;;, with the appropriate substitution for m we get: 


(8) 


40) 


Wh 


1— cos lop T 


Using Eq. (8) we can determine the number of turns in the load winding required for any given value of top, 


The preparatory pulses may be applied to the magnetic element by means of various remote-control de- 
vices. Figure 4 shows the circuit of the element I, the source of the preparatory pulses. The output winding of 
this element is represented by Wy, Rpy and Rpg represent the dynamic resistances of rectifiers 1 and 2, and R, re- 
presents the resistance present in the preparatory circuit, its basic function is to limit the return information [4]. 


I yi 
“,3 W, > 
My, Rr 
pulses 
RifRey The diagram given in Fig. 4 is the outline 
9 — of a circuit which serves to transfer the current from 
Umasinwt Y, 24 core I, which can be considered the sending ele- 
re ment, to core II which functions as the receiving 


element. This type of circuit is well known and 
has been described in the literature [4]. 


It has been shown in [4] that for cores I and 
II of the same size and with identical characteris- 
tics, the operation of the sending core must occur more slowly than the preparation of the receiving core, i.e., 
it is necessary to satisfy the conditions to, ; > tprep Ilr where tprep 11 represents the time required to prepare ele- 
ment IL. 


Fig. 3. Circuit of pulsed magnetic element 
(without feedback paths). 


As a preliminary step, let us assume that t,1 = threp]- Im this case, during the preparation of element II 
of the preparatory circuit, we have 


t 
= Tua + Pian rdt + (9) 


Using the piecewise-linear approximation for the volt-ampere characteristic of rectifier 1, we get 


UB, = Uara + 


where Ug, ; is the threshold voltage in the forward direction of rectifier 1, Rpg is the dynamic resistance in the 
forward direction of rectifier 1, and i, is the current resulting from the electromotive force ey. 


Substituting the value of up, in (9) and integrating, we get 


1 (Ry +R 


Applying the known condition Wy = 2W;, and solving the latter equation for Wr, we get 


Marr fon 10° S47 opr Al + Ry top 10° 


During the operation of element II the preparatory winding W; may send the winding W), return informa~- 
tion which can partially prepare the sending element I, and in this mamner lead to the presence of conditions 
which produce interference. 


Let us assume that there is a tenfold weakening of the voltage applied to the Wy, winding when the second 
element operates. This condition is approximated by the relations in i,R, = 10ugy = 10ugrs + 10ipRpg, where ir 
represents the current flowing through resistance R;. Neglecting the term 10ugrg, we obtain Ry = 10Rpg. Since 
Wx = 2W,, then the ratio of the change in current due to partial preparation of the core as a result of the return 
of information to the change in current resulting from complete preparation is approximately equal to 1/20. Af- 
ter the substitution of Ry = 10Rp, in (9a) we get 


Warr 10° 7 tary 10 \2 Ae ly (Ruy + 10 10° 


From this expression we can determine, for given types of rectifiers 1 and 2, the magnitude of W, and the 
resulting Wy. These calculations give us the limiting case top] = tprep- However, for this case the operation 
of the circuit shown in Fig. 4 is not stable, since complete preparation of element II is not assured should the 
supply voltage be lowered. In order to insure stable operation of the circuit we will decrease the magnitude of 
Ry, so that 


Ry = (5-7) (10a) 


Then the ratio of the complete to the partial preparation is approximately equal to (5-7) W/W, = 10-14, 
which is perfectly acceptable. This assures the condition topI > tprep II: We wish to note that if the parameters 
of the preparatory circuit (Rr, W;) are calculated on the basis of operation under the condition top; = 0.5T, then 
the circuit operation will be stable if to, > 0.5T. 


As a matter of fact, it follows from Eq. (1) that 


‘opi 
e,dt = S = const, 
0 


where S is the area bounded by the curve e; and the x-axis, as shown in Fig. 5. Therefore, we may write, that 
for top = 0.5T 


S, + Sy = +S, =5, 


where the areas S, and S*; are proportional to the change in inductance of the element IJ, in the preparatory pro- 
cess of reversing the magnetic polarization of element I, during the time top1 = 0.5T, where top1< 0.5T, and 
the areas S, and S}, are proportional to the energy absorbed by the resistor Rr during the time top 7 = 0.5T, where 
top < 0.5T. 


It is clear from Fig. 5 that S$ < S, and therefore S'; > S;. Therefore, we see that if the area S, is to be of 
sufficient size to insure the complete preparation of element II, then the area S'; will be more than adequate to 
insure the complete preparation of this element. 


The current in the input winding W, (Fig. 3) is a maximum when the element is in the “nonprepared” state, 
i.e., when the preparatory pulse is absent and the inductive reactance of winding W, is practically zero. For this 
case, the rectified current in the input circuit is given by 


where Rp is the additional resistance in the input circuit and U 
is the effective value of the sinusoidal voltage applied to the in- 
put circuit. 


Given the diameter of the conductors of the input winding 
Wy, we can determine the value of the current I) so that we will 


not exceed it. 
It is obvious that the condition 
Fig. 5. Curves of the voltage in the 1, =X <Ip (11) 
preparatory circuit during the time in- D 
terval top = 0.5T, where top < 0.5T. ae 
Using the above condition we can determine the additional 
resistance 
2 
(12) 
For the input circuit 


The current i, is determined by the expression 


i; = ig — 


Ww 
— wt (14) 
Substituting the value of i, from (14) into (13), and integrating both parts of the resulting expression be~ 
tween the limits 0 and to, (operating time of the element), we get: 
t 


+ Ryidt — — Ryniedt. (25) 


Performing the integration, we get 


2U V2UH, t,t 


W,(Ry + 17-40-47" W,) 


BLS 4 1.7-10-"W2 (16) 
? 


+ 1,7-40-W,, 


As a result of Eq. (16) we will consider the magnetizing current constant and equal to Hol, /0.49W, for 
the time t = 0 to t = top. 


Equation (16) may be placed in the form 


BU (1 — top) = kW +h, + (17) 


where the coefficients k';, k, k's, k4 and k's are constants. 
Using Eq. (17) we can determine W, given U and top, or given W, and top we can determine U. 


Calculation of the parameters of a two-core pulse relay with resistive load. The circuit of such a relay is 
shown in Fig. 6. 


Fig. 6. Circuit of a two-core (RHL) pulse relay. 


The method derived above may, with a few corrections, be used for the calculation of the given relay. 
The number of tums in the load winding of each core is given by the expression 


U,, 
W,= 


6,28B,S, 


We see that in this case the expression in parentheses, under the square-root sign, is multiplied by the factor 
2/T, rather than by the factor 1/T,as was the case for the element without return paths. This is due to the fact 
that current flows in the load winding during both halves of the input voltage cycle. 


The number of turns in the start winding W; and the resistance R; are determined, as before, from Eqs. (10) 
and (10a). 


In order to determine the number of turns in the intermediate connecting windings Woo, and W4on, we 
may use Eq. (9a) with Roop substituted for R,. 

By analogy with the case described above, it is necessary, in order to assure stable operation of the relay, 
to choose a value for the resistance R',,,, that is somewhat lower than that used in calculating Woon and Won! 
R'con = (0.5-0.7) Reon: 

In order to determine the input parameters Rp, U and W, of the relay we may use Eq, (12), and an expres- 
sion analogous to (17) may be used to take into account the supplementary effects due to the ampere-turns of the 
intermediate connecting winding. 
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Analysis of the Operation of a Single-Core Pulse Relay with Resistive Load 


Let us analyze the operation of a relay with a resistive load, the circuit of the relay being represented by 
Fig. 7, where We, R, and C represent elements of the feedback path. 


In order to simplify the analysis of the relay operation we will neglect the forward resistance of the recti- 
fiers and the winding resistance. In addition to this, we will also neglect the influence of the start and stop cir- 
cuits. These assumptions are permissible since the influence of the in- 
dicated factors upon the operation of the relay is negligible. 


tS During the operating half-period 
Ww, 
where 
du, 
W, pel i, = tp + ic, ic =C—-, 
Ww, 
Up axsinwt U, 
re We have found, by means of experiments, that the most stable 
operation occurs in the case where the relay operates during the time 
top = 0.25T, wt = 90°. 


Then during the time that the relay is in the operated condition, 
€, = Ug and therefore, 


W, W, de. W,, ¢ 


The following relation holds for the relay input circuit 
u = igh + Uy ‘ (20) 
Substituting the value of iy obtained from Eq. (19) into Eq. (20) and making use of the relations 


w 
we obtain the differential equation 


Solving this differential equation for the initial conditions uy = 0 and t = 0, we get 


U 
u= Sin ot 008 wt) + 


KC kip | CKp 
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After some manipulation Eq. (22) may be placed in the form 
U, 


= sin (wt —aretg ©) + 


We see from Eq. (23) that the voltage in the relay winding consists of three parts, the first part is a sinus- 
oidal voltage of the same frequency as the input voltage but with a different phase and amplitude. The second 
part contains the factor e~P* and represents an exponential voltage. The third part represents a constant voltage 
and is the result of the magnetizing current in the core. 


During the operating period of the relay (Fig. 8) 


ty 
| ujdt = 2 4W,-10"*. 


If we substitute the value of u; from Eq. (22) into the above equation the equation may, after some manipu- 
lation, be placed in the form 


gx? U.__wsin et; 
= 2B (24 


At the end of the magnetization (operation) the voltage across the capacitor increases to the value uc; 
(Fig. 8). After this the current changes in the core and consequently the capacitor begins to discharge through 
the winding Wy. 


The discharge current tries to prepare the relay for its next operation but, since current is still flowing in 
the input winding W, up to the time t = 0.5T, and the ampere-turns iyW, > ic dischW2, this preparation does not 
occur. The inductive reactance of the winding W, is zero for the time t — ty, and the discharge time constant 
is RgC. During this time, the capacitor voltage decreases, following the rule 


uc=ene (25) 
For t = ty 


4R,C 


UCg = & 


(26) 
where e,, represents the electromotive force in the feedback winding during the capacitor discharge for t = ty = top 


At the time t = t, the current i; ceases, and the winding “becomes part of” the discharge circuit, from this 
time on and until the polarity of the core is completely reversed, the discharge current will be constant. 


In order to insure stable operation of the relay during the changes in the input voltage (decreases), provision 
must be made earlier to insure an adequate voltage supply. With this aim in mind, we require that the core re~ 
turn to its initial state during the time t, —t, = 1/sT, where ts < & = T. 


From Fig. 8a we get 


a b 


NS § 


Fig. 8. Voltage curves for the feedback circuit. 


where imag represents the magnetizing current, which is equal to the capacitor discharge current during the re- 
turn of the capacitor to its initial state. 


AUc= = inagl!s — 


Since we have assumed that t, —t, = ¥5T, then AUc = imagT/3C. In order that the element be fully pre- 
pared during the time t, — ty, it is necessary that 


U 
(te — ts). (27) 
The latter expression becomes clearer if we refer to Fig. 8b, where the area S represents the number of 
volt-seconds required for the complete reversal of the magnetic polarization of the core. 
Using the basic Eq. (1) we get 


ts 


| = 2B,S 


After substituting the values of Uky and Ujg and rearranging the resulting expression, we get 


from which 


2B,S + 4Rsimagts — ta) + ‘ae 


eg, = 


t3 — ty 


If, in the latter equation, we let ts = T, then we get the minimum value of the electromotive force in the 
blocking winding that is necessary for the complete preparation of the core. 


If, in Eq. (28), we make the substitution shown below 
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T T WwW 
en = —Un 


where u,; is the voltage in the input winding at the end of the reversal of the magnetic polarity, this may be ob- 
tained from Eq. (22) if we let t = t;, then we get 


iy 


At the time t, the core is saturated. Therefore, during the time interval \ — t, the inductance of the W, 
winding will be zero. As a result, the capacitor will continue to discharge with time constant 4RgC. The voltage 
across the capacitor during the discharge time t — t, is given by the equation 


then at time t = &, we have 


Uca= . (31) 


After the time t = t,, the second input pulse appears and the relay operates, Fig. 8a. While this is occurring, 
and up to the time t = tg, there is joint action of the input and blocking circuits, this action may be described by 
means of a differential equation of the first order, the solution of the equation consists of three terms as shown in 
Eq. (23). At the time t = ts, the capacitor starts to charge, after which the process repeats itself, according to 
the description given above. 


The expressions (23), (24), (29) and (31) may be used to determine a series of relay parameters. 
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MACHINE ANALYSIS OF SWITCHING CIRCUITS 


P.P. Parkhomenko 


(Moscow) 


Algorithms for the mechanization of the structural analysis of switching circuits 
and for the sequential determination of the action of their elements are presented. The 
realization of the algorithms in a special logical machine is described. The machine 
performs a complete analysis of switching circuits of practically any complexity much 
better and mare rapidly than man does. The results of the analysis are recorded man- 
ually or automatically. The machine may be used for the analysis of switching circuits 
with solid-state elements and for the solution of logical and other types of problems. 
The principles involved can be used for automatic quality control in the production of 
relay devices and other electrical apparatus. 


1. Introduction 


Machines composed of relays are widely used in all branches of technology and national economy, they 
are often extremely complex and a great number of qualified intellectual workers are occupied in the design of 
these machines and their maintenance after assembly. The mechanization of the design of relay machines would 
free these workers from many monotonous and tiresome operations, would sharply raise the productivity and qual- 
ity of their work and would provide them with the possibility of investigating newer and more difficult problems 
of relay technology, the solutions for which, in spite of the great success and stormy development of the theory 
of relay machines, are not at the present time within their ability. Also, as an end result, the specialization of 
machines for the solution of the problems involved in the synthesis and analysis of relay devices would greatly 
influence the further development of science and technology. 


Efforts have been made both here [1, 2, 3] and abroad [4, 5] to deyise a machine that will facilitate and 
accelerate the process of switching circuit design. Papers [3,5] which are devoted to the problem of the mech- 
anization of the synthesizing process will not be examined here, 


The switching -circuit analyzer of Shannon and Moore [4], has been designed by the authors for use in con- 
nection with the synthesis of two terminal switching circuits made up of the contacts of at most four relays. This 
analyzer allows us to determine whether the circuit being investigated satisfies the specified conditions of opera- 
tion, it is simple and consists of a small amount of apparatus. The principles upon which the analyzer is based 
are not limited to two-terminal networks or to four relays, however, for example, each addition of one relay to 
the circuits being analyzed would double the size of the machine and the quantity of apparatus, The analyzer 


cannot be applied to relay circuits involving many contacts or to circuits involving the relay windings as well as 
the contacts. 


Tsukanoy's matrix analyzer [1] mechanizes part of the operation by matrix analysis of the switching cir- 
cuits. The time required for the analysis depends upon the number of relay contact chains in the circuit being 
analyzed rather than upon the number of relays in the circuit as in the Shamnon~Moore analyzer. 


*This machine has been developed by the author under the direction of Professor M.A. Gavrilova. 
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The results obtained by use of the analyzer cannot be applied directly, without supplementary processing, 
to the design of a circuit in the general forme that fulfills the desired conditions. The investigation of multi- 
contact circuits by means of the matrix analyzer is inconvenient and also it is not sufficiently mechanized. 
Complete automation of the process of analysis is impossible because of the method that has been chosen to in- 
troduce the circuit under analysis into the machine (the structural formulas for the circuit chains are written in 
chalk, between the loops on the board located in the front of the analyzer. All of the three versions of the mat- 
rix analyzer that have been constructed (the one using electromagnetic relays, the one using step-by-step switches 
and the one using perforated cards) do not permit the analysis of a circuit involving more than six loops, and in 
no way differ from one another in the functions that they fulfill. 


Fig. 1. 


Rodin‘s electronic analyzer [2] is an electronic analog of the Shannon and Moore analyzer. The electronic 
analyzer has the capacity required for the analysis of a six-relay circuit. This analyzer is also the first fast-acting 
machine of its kind. While it is possible to increase the capacity of the analyzer, this is very difficult due to the 
resulting large increase in the size of the machine. 


In contrast to the existing machines, the machine that is described in this article, the IAT AN SSSR, does 
not have the principal limitation, that of capacity, in addition it can work with both single-cycle or multiple- 
cycle circuits, and also with switching circuits of the general form that have already been assembled and built. 
The machine verifies the circuit operation, either for all combinatorial states of the circuit elements, or for cer- 
tain definitely designated circuit combinations, for multiple-cycle circuits the machine indicates the sequential 
action of the circuit elements directly, and the results may either be recorded manually or they may be printed 
automatically by the machine. The machine may also be used to solve a number of other problems of a tech- 
nical-logical character. The exterior view of one of the first models of the machine, together with the printing 
apparatus, is shown in Fig. 1. 


A machine with the capacity for analyzing twenty element switching circuits consists of 180 relays RKN, 
six step-by-step switches of the type SHI-25/8,** 500 tri-positional switches and 120 signal lamps. The overall 
size of the machine is 1800 x 1200 x 350 mm. If we desire to increase the capacity of the machine so as to en- 
able it to analyze the operation of switching circuits containing one more element, then it would be necessary 


*By the general form of a circuit we mean a form in which the relay contacts (or other reacting organs of the 
circuit elements) are shown as part of either a shunt or series connection. 
** Instead of the indicated step-by-step switches three reversible switches of the type RSHI-25/8 may be used. 
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to add eight more relays, two more mounting boards of switches and a proportionate amount of other apparatus. 


We wish to note that it would simply not be realistic to build an analyzer of the type described in [4], that 
would have the capacity to analyze twenty-element switching circuits, since this would require switches with 
1,048,576 contacts just to account for the states of the constituent elements, not to mention the rest of the ap- 
paratus that would be required. 


Reference [6] gives a short description of the first version of the machine IAT AN SSSR that was constructed 
and its operating conditions. At the present time a more perfected type of machine, built in a more economical ' 
block version, has been developed and it is this machine that will be considered in the following discussion. 


2. Mechanization of the Analysis of Circuit Structure 


In the general case, switching circuits consist of receiving, intermediate and executive elements, each of . 
which has reacting and performing organs. The problem involved in analyzing the structure of a switching cir- 
cuit consists in determining the operating conditions of the reacting organs of its intermediate and executive ele- 
ments and their dependence upon the state of the performing organs of the receiving and intermediate elements. 


Let us designate the reacting organs of the various elements by capital letters, and their performing organs i . 
by the corresponding small letters. Let us use for the performing organs p of the receiving elements of the circuit, 


the designations 4;,..., aj,..., @p, where a is the same as aj or aj, let us use for the reacting organs q of 
the intermediate elements of the circuit, the designations X,, ... , Xj, . . . , Xq and for the performing organs 
of the intermediate elements of the circuit, the designations Xs... , Xj,...» Xqs and let us use for the react- 
ing organs r of the executive elements of the circuit, the designations Z;,..., Z,..., Zg. In order to be 


definite, we wish to state that this paper will consider basically switching circuits that are constructed of two- 
position elements. However, taking into consideration the generality of application of relay circuits, the results 
obtained may be,to a great extent, applied to circuits composed of “yes-no” elements of any type. 


When classified according to their principles of operation, switching circuits may be single-cycled (circuits 
of class O), in this case q = 0, or many-cycled (circuits of class M). From the structural point of view the cir- 
cuits may be classified as parallel-series (circuits of class P) or circuits with bridged connections (circuits of 
class N). 


In order to solve the problem of the mechanization of the structural analysis of switching circuits, it is 
necessary first of all that we determine a method of analysis that can be applied to circuits of all classes in the 
general form. The degree of formalization of the chosen method must permit the construction of an algorithm 
which permits a sufficiently simple realization of the problem in the machine, An investigation of the existing 
analytic, graphical, matrix and other methods of analysis shows that the most suitable method of analysis for our 
purpose is the method of resolving the circuit into its constituent units [7]. ; 


In the proposed form [8] this method is applied only to the structural analysis of two-terminal contact net- a 
works of Class II (p > 0, q= 0, r= 1). The method is based on the fact that any Boolean function p may be re- 
presented as the sum of terms, each of which consists of the product of all the variables p or their negatives, ; 
taken in various combinations; } 


(@4,..., apy) = hi Il (1) 
i=! 
Here the aj are the variables, the hj are the coefficients of the variables, the coefficients assume the values 
0 or 1, £ is the summat.on sign and II represents the product (multiplication). 


Analysis of the conditions required for the operation of the two-terminal relay circuit by means of (1) leads . 
to the equation, which sets forth the combinatorial state* of the receiving elements, 


Pp 
i=1 


*We will designate the combinatorial state of the circuit elements by means of the product of symbols as in (1) 
or (2). 
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for which a closed path exists between the terminals. These combinations correspond to hy = 1 and represent the 
complete elementary chain of the two-terminal contact circuit. An analogous set to that of Eq. (2), corresponding 
to hy = 0, gives the results of the analysis for nonoperation. 


Extension of the analysis (1) to two-terminal relay contact networks of class N is based upon the fact that 
any such two-terminal network, when analyzed according to the initiating or terminating elements [7], may be 
converted to an equivalent two-terminal network of class P. 


For given conditions of operation the elements of a circuit, which is of the general form, usually satisfy 
the formula [9] 


Fr = (3) 


where the Z; represent the reacting organs of the elements being analyzed, f[Z},] represents the chain switched 
in as a result of the action of the Zy, and f{Z,} is the inverse chain of parallel Zy,. Since the chains f[Z,] and 


f {Z},} represent two-terminal networks, then instead of (3) we may write 


quire that we consider the set of elements yy, for which 


{ha} {hi} = 4 (5) 
and the conditions that must be fulfilled for the nonoperation of the element Z,, require that 


(hi) {hi} = 0. (6) 


The physical analysis of switching circuits by means of the division of the switching circuit into its con- 
stituent units depends upon the determination of the state of the relay windings and the state of the relay contacts. 
For this it is essential that changes in the state of the contacts be entirely independent of the state of the relay 
windings. in all cases for of elas O, but is physically contradic- 
tory for circuits of class M. 


wwe dnslyme the cleouits of clam M, we see that it is postible to subdivide them into their constiment tilts 
only if during the analysis we turn our attention away from the existing dependence between the relay windings 
and the contacts of the intermediate relays. In other words, in the course of circuit analysis we cut the return 
paths that exist, for this is required to change the circuit to one of class O. Let us call the latter circuit, the 
single-cycle circuit equivalent of the original class M circuit (or, more briefly, the equivalent circuit). It is 
clear that the receiving elements of the equivalent circuit G1, . . . , &p+q are the receiving elements ay ..., 
G, and the performing organs X, . . . , Xq of the original circuit of class M; the performing elements of the 
equivalent circuit ¥,,..., ¥q,, serve as the reacting organs of the intermediate elements X;,.. . , Xq and the 
performing elements Z,, .. . , Zy of the original circuit. We will call the receiving elements of the equivalent 
circuit the input elements and the performing elements of the equivalent circuit the output elements. 


For the single-cycle equivalent circuit, Eq. (4) holds, provided the proper symbols and indices are sub- 
stituted: 


=0 i=1 i=1 
where k=1,...,4+14. The expansion of (7) holds for the most general case, for q = 0 we get (4), if the shumt- 


ing chains are absent we get (1); p = 0 leads to the autonomous circuit (10), while r = 0 leads to a circuit without 
performing elements. 
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The mechanization of the analytic process assumes that the 
circuits to be analyzed are introduced into the machine in some 
manner. The method of introducing the circuit into the machine 
influences, to a significant degree the complexity of the algorithms 
realized by the machine, The method chosen for use with the ma- 
chine described in this article, the introduction of a physical model 
of the circuit to be analyzed, is the most direct and natural method 
of introducing the required data, not only does it permit us to realize 
in a sufficiently simple manner the algorithm of the structural 
analysis and, as developed below, also an algorithm that enables us 
to determine the sequential action of the elements in a switching 
circuit of class M, but it also insures the undisputed pre-eminence 
of specialized machines, even those constructed of electro-mech- 
anical apparatus, over the universal rapid-action computers for the i 
investigation and analysis of switching circuits.¢ 


The presentation of circuits of class M in the form of single- 
cycle equivalents is necessary in principle in order to model the 
circuit physically, it insures a uniform process of structural analysis 
that is independent of the class of circuit being analyzed, and it is 
also useful in connection with the investigation of questions concern- 
ing the determination of the sequential action of the circuit ele- 
ments. More than this, an understanding of the operation of the 
single-cycle equivalent circuit permits us to apply all the methods 


of analysis and synthesis developed for circuits of class O to circuits , 
of class M. 


Taking into account the fact that the circuit of class O is, in 


the sense determined above, a particular case of the single-cycle equivalent circuit, we will, in the following 
analysis, speak in terms of the latter circuit. 


The physical model of the circuit may differ to a significant degree from the relay circuit intended for in- 
dustrial use in its construction, in the parameters of the apparatus used, the time scale, and other such factors. 
However the model must preserve all the circuit characteristics essential in any given case for the aims of the 
analysis. This, first of all, refers to the circuit structure, i.e., the relationships between the elements and the 
character of the electrical connections between them, and in many cases to the mutual parametric dependence . 
that exists in circuits of the general form. The basic elements used in the modeling of switching circuits consist 
of relay contacts (make and break contacts) which take the part of input elements, the relay windings which are 
used in the capacity of output elements, and the conductors which provide the necessary electrical connections 
in accordance with the structure of the equivalent circuit. In addition to this, it is advisable that the model con- 
tain among its elements resistors and rectifiers. Of course the model must also possess a power supply. ; 


Due to the lack of space we will not discuss such general questions as how to determine the capacity of 
the model, how to determine suitable parameters for the elements of the model, the design of the power supplies, 
the electrical connections, and other such details that may be involved in the modeling of switching circuits. We 
will limit ourselves here to a short description of the decisions that were made in the construction of the machine 
that is being described in this paper. The input elements of the circuits to be analyzed are represented in the 
model by the contacts of the receiving relays Cj, and the output elements are represented in the model by the 
windings of the performing relays Kj. The contacts and relay windings are connected to sockets which appear on . 
the plug-in board of the machine. Two sockets are provided for each electrical element of the circuit being 
analyzed to facilitate the multipling of connections associated with each loop of the circuit. Each relay C; is 
provided with three transfer, one break and two make contacts, and each relay Kj has three identical, separate 
windings, The machine is provided with an equal number of receiving and performing relays, the basic factor 
in determining its capacity. In this manner the windings of the Cj relays may control the action of the input 


elements and the reactions of the circuit being analyzed may be placed under the control of the contacts of the 
Kj relays. 


*A convincing comparison in this respect is given in [4]. | 
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The windings of several auxiliary relays Dy, each with three windings, are also brought out to the plug-in 
board. These relays are used as supplementary relays, providing additional contacts should the number of C; re- 
lay contacts prove insufficient to model a specific surface undergoing analysis, or they are used for other special 
purposes. Similarly, auxiliary resistors, rectifiers, etc., are also provided. Direct current is supplied to the ma- 
chine by means of several supply sockets that are provided both in the upper and lower parts of the plug-in board, 


p, 


Fig. 3. Block diagram of the machine used for struc- 
tural circuit analysis. 


Electrical connections required in the model are set up manually on the plug-in board of the machine by 
the use of double-ended connecting cords, with a jack at each end, the jacks being plugged into sockets as re- 
quired, in order to construct a model of the circuit being analyzed. An example of a model circuit, set up on 
the plug-in board of the machine is shown in Fig. 2. 
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Fig. 4. 


The realization in the machine of the algorithm of the structural analysis of the circuit is accomplished 
in the following manner. A special device, the constituent generator GC (Fig. 3), is used to feed all pos- 
sible input combinations, in time sequence, to the receiving relay RR shown in the block diagram of model 


M (consequently also designating input elements to the equivalent circuit. When a circuit is being 
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analyzed to determine the required operating (nonoperating) conditions, should a certain output combina - 
tion resulting from the action of the receiving relays of the model lead to the operation (nonoperation) of even 
one of the performing relays, PR, then the machine will stop and present the particular result of the analysis, 
which includes the input combination {ay . . . , &p4q} and the corresponding output combination {¥,..., 
Yqsr} of circuit elements. If the results are to be copied down manually they are obtained by means of lamp 
signals, after which the operator allows the machine to proceed. If the printing device is used, the machine pro- 
ceeds with the problem automatically after it has finished printing the partial result. 


The complete results of the analysis are presented in a table, the lines of which consist of a pair of com- 
binations; 


Here m =1,..., p denotes the number of the result in the sequence and p denotes the number of partial 
results of the analysis. 


The method of analyzing a switching circuit by separating it into its constituent units possesses many ad- 
vantages, To use the results of the analysis we simply compare them with the results obtained by synthesis as 
shown in truth tables for circuits of class O, or in action tables for circuits of class M. This type of analysis al- 
ways results in a definite answer as to the relative merits of two circuits under comparison; allows us to deter- 
mine the influence, upon the state of several elements, of the action of very fast, short-time duration chains 
occurring due to a change in one of the cycles of a class M circuit; simplifies the uncovering of previously un- 
foreseen work cycles that may occur during circuit operation, and so on. 


However, this method possesses a major disadvantage, in that the length of the analytic process is doubled 
for each additional input element that is added to the circuit undergoing analysis. This defect cannot be re- 
moved by constructing the machine entirely of elements that do not possess inertia because if we do not use a 
physical model of the circuit in our analysis, then a sharp increase occurs in the number of operations required 
to determine the number of closed paths that may exist in the circuit, not to mention the the difficulties in- 
volved in this case, in coding the circuit for presentation to the machine. However, this defect can be compen- 
sated for to a significant degree by building the circuit of the constituent generator in such a manner that n, of 
the input elements are kept constant for a series of inputs during the proce:s of analyzing the input combinations 
In this way the time interval that is required for the introduction of the various possible combinations is shortened 
21 times. We have done this in our machine by the use of the three-position switches XY, which are shown in 
the circuit diagram of the constituent generator (Fig. 4). With the switch in the upper position (x contacts closed) 
the paths which control the action of the corresponding receiving relays of the block diagram M are closed, 

With the switch in the lower position the circuit connections which control the receiving relays are open. Con- 
tacts u and u belong to the impulse generator relay; they control the operation of the constituent generator. The 
design of the constituent generator circuit is obtained analytically; this circuit insures a sequential presentation 
of the constituents in which each succeeding combination differs from the preceding combination by a change in 
only one factor. 


3. Mechanization of the Sequential Operation of the Circuit Elements 


The determination of the sequential action of the circuit elements is the most important stage in the analy- 
sis of the operation of class M switching circuits. 


Now, the operation of the switching circuit is, in the end, determined by its structure.A complete under- 
standing of the structure may be obtained by an analysis that is based upon the subdivision of the single-cycle 
equivalent circuit into its constituent units, and the results of this analysis may then be used to determine the 
sequential operation of the circuit elements. The solution of the problem consists in the construction of an algo- 
rithm that will enable us to transform the results obtained from the structural analysis (8), as Listed in the form 
of a table of combinations, into a table that will give the required sequential information. 


It is not difficult to show that such an algorithm may be built upon the basis of the obvious fact that for 
the single-cycle equivalent circuit the following relations hold; 


*This follows from the fact that the number of constituent units is equal to 2", where n is the number of input 
elements in the circuit. 
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where t and t + r are numbers designating succeeding steps in the resulting sequence table. Therefore, given 
the input combinations initially, and knowing the order of operation of the inputs upon the receiving elements, 
it is possible, with the aid of the algorithm,to uniquely determine the possible sequential operation of a class M 
switching circuit. Naturally, the sequence table that is obtained by this process is realizable.* 


In spite of the high degree of formalization connected with the indicated algorithm, its realization is con- 
nected with significant difficulties due mainly to the fact that a memory with an enormous capacity is required 
to store the table (8). We were forced by the latter circumstance to seek a more rational method of solving the 
problem. 

Let us assume that the circuit of class M that is to be analyzed is represented by a physical model of the 
corresponding one-cycle equivalent circuit. Let usalso assume that the initial combination of input elements,as 
well as the program of external influences upon the receiving elements of the circuits, have also been provided. 
In order to represent the principal logical operation of the circuit in the machine by means of an algorithm, 
which is to be used to determine the sequential action of the circuit elements, we will introduce the following 
designations: M to represent the physical model of the single-cycle equivalent circuit, B to designate the circuit 
block of the machine that represents the initial combination of input elements and the changes in the settings of 
the receiving elements in the steady state of the tircuit being analyzed, O to designate the block representing 
the return paths, which insure the transfer of the settings to the performing relays of the model, the latter relays 
functioning as the reacting organs of the intermediate elements of the circuit, corresponding to Eq. (11), D to 
represent the block that provides a time delay to take into account the operating or release times of the circuit 
elements, "Out" to represent the output of the machine which gives the settings of the circuit elements as time 
steps in a sequence table, Y, (t) to denote information obtained from block A at time t, G(Y, (t) - C) to repre- 
sent the transfer of the information Y, (t) to block C, C[¥, (t) ~ Y¢(t)] to represent the operation of block C 
that results from the information Y, /t), changing the state of the letter to Y¢(t), p(A,C) to represent the logical 
conditions that determine the block (A or C) to which it is necessary to forward instructions, w to represent the 
logical conditions that are always false. 


Using the accepted conditions (11), we obtain the logical scheme for the algorithm in the following form: 


| BG (p(t) out M)M (Ya(t)~Yu(t)) | G 3) 
(Yp(t)—> out, 0)O(Yp (t) ~ Yo(t)) 
G(Yo(t)—> M, B) p(M, B) t (22) 
The block diagram of the machine, shown in Fig. 5, is based upon the principles developed in the logical 
analysis (12). The machine operates in accordance with the principles outlined as follows. 
The first combination of input elements for the circuit being analyzed 


is placed into the machine with the three-positional switches** of the input block B in the upper position, and 


*The table of sequential actions for class M circuits is realizable, if and only if, it does not contain at least two 
steps tandt+ At, At > 1, with the same combinations of receiving and intermediate elements, each of which is 
immediately followed by steps requiring different combinations of the intermediate and performing elements. 
These are the same three-positional switches XY that are shown in the block diagram of the constituent gen- 


erator (Fig. 4). 
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transferred at the same time both to the receiving relay RR of the physical circuit model M and to the output of 
the machine, where it is placed in the zero time position in the output table of results. The input combination 
(13) leads to the output combination 


which is obtained from the performing relay PR of the physical circuit model M, this output combination then 
enters the time-delay block D, The latter operates upon the combination (14) in such a manner as to prevent 
any change in the elements of the above combination (prevents the operation or release of any of the delayed 
elements) during the time interval sr of the delay. After leaving the time-delay block D the output combina - 
tion together with the combination 


{ay,.. (15) 


of the remaining receiving elements, which have not undergone any change, is fed into the output of the machine 
as the first line in the output table. At the same time, the combination of intermediate elements 


(16) 


enters the feedback block O which feeds the combination to the corresponding receiving relays of the model M. 
In this manner, if the circuit has not reached a stable condition, its receiving elements are fed a new combination 
that differs from the combination (13) and the process repeats itself. 


M 
8 A, Aa} eve 


RR PR At é sr 


4, Ag 4tésr 


0 


Output 


Fig. 5. Block diagram of a machine used to determine the sequential opera- 
tion of the elements of a circuit. 


The stable condition of the circuit is characterized by the equations 
{X,, Xahi4+e {X,, Xahts (17) 


{Z,, Zr}t4+ {Z,, (18) 


whose action is controlled by a special loop of the machine (not shown in the block diagram of Fig. 5) and cor- 
responds to an untrue result of the logical conditions p(M,B) of (12). When conditions (17) and (18) occur, the 
machine stops and, in order that the process continue, it is necessary in the general case, to change the setting 
of the receiving elements of the circuit. This is accomplished when the operator sets the switch of input block B 
to the lower position. 


The moment during which the machine completes the process of obtaining the table of output values is not 
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automatically determined because this would require that the machine be provided with a memory of unjustifi- 
ably large capacity. Practically speaking, the operator can easily determine this moment since, as a rule, it 
coincides with the stable state of the circuit which occurs at the same time that the machine stops. 


In order to facilitate manual recording of the results, the table of output values is placed, in parts, on an 
illuminated board. The machine may also be constructed, by means of an additional attachment, so that it 
prints the steps of the resulting table of combinations automatically. 


4. Other Applications of the Machine 


Analysis of switching circuits constructed of solid-state elements. The ‘possibility of using a machine that 
has been constructed to analyze relay switching circuits for the analysis of circuits constructed of solid-state 
elements is based upon the fact that any such circuit may be represented by an equivalent relay circuit, which 
may be placed into the machine and investigated in the usual way. 


The machine is constructed in such a manner that it is possible, for example in an electron-tube circuit, 
to determine the operating conditions at any point of the circuit for all possible combinations of the input con- 
ditions, the operation of the machine being limited only by the action of control voltages applied at the outputs 
of the circuit. 


Comparison of relay circuits. The two circuits that are to be compared are placed into the machine in 
such a manner that they make use of the same performing relays of the machine; however, a separate winding 
is used for each circuit, the windings being placed on each relay,series-aiding. Therefore, during the structural 
analysis, the machine will stop only for those constituent units of the circuits being analyzed which are present 
in one circuit and absent in the other. 


The process of comparing the two circuits permits us to examine each circuit in detail, and enables us to 
simplify the circuits, not only by shortening the circuit paths by the omission of individual relay contacts, as is 
done by the analyzer of Shannon and Moore, but also by the use of other methods, such as the combination of 
similar elements, and so on. 


Solution of logical problems. The machine may be used to solve many logical problems, since any such 
‘problem can be presented in the form of conditions, in which the variables and their negatives are bound to- 
gether by means of the relations “and” or “or.” These relations correspond either to a series or to a parallel con- 
nection of a relay circuit. For this reason a logical problem may be presented in the form of a relay-contact 
circuit, placed in the machine, and solved in the usual manner, i.e., the action of the circuit may be investi- 
gated for all possible input combinations of its variables. 


In principle, the machine may be used to solve all problems, which can be presentec either directly or in- 
directly in a form suitable for placement into the machine. 


Automatic quality control in the production of relay devices. The machine described in this paper may 
be used to automatically detect errors in the assembly of various relay blocks, or other irreparable circuit errors 
occurring during the manufacturing process. We wish to point out that the first models of the machine were not 
designed with this purpose in mind, and for this reason they possess a series of limitations in this respect (such as 
the presence of only one source of supply, the absence of automatic means of entering the operating conditions 
of the circuits under test into the machine, etc.). However, the principles upon which the operation of the ma- 
chine is based permit it to be used successfully to control the quality of series-produced electrical apparatus and 
devices, if various specialized trial and test stations are established in the production line. 


SUMMARY 


Experience gained during the use of the machine described in this article for the solution of the practical 
problems involved in the analysis and design of relay-switching circuits confirms its effectiveness. We have 
found that the use of the machine is a great advantage in the analysis of multi-element, multi-step circuits. 


The machine described in this article, which is used to analyze the operation of relay circuits, received 
high praise in the foreign technical press, and from among the other automatic devices it was chosen for the 
Grand Prize at the exposition held in Brussels in 1958. 


At the present time we find ourselves at the stage of creating universal machines that are capable of de- 
signing relay circuits. Machine-circuit synthesis, a process of extremely great value, will undoubtedly require 
machine analysis as the final step of circuit design. In this case, the machine that is used for the circuit analysis 
may be substantially simplified if the data that represent the circuit to be analyzed are automatically introduced 
into the machine, since this will permit us to remove the circuit loops that are required in connection with the 
manual control of the machine occurring during the process of obtaining and verifying the output table. The 


most likely way of introducing the circuit to be analyzed into the machine seems to be through the use of punched 
cards. 


The author wishes to thank Prof. M.A. Gavrilova for the valuable observations and advice which she gave 
him while he was working on the design of the machine, which was to be used for the analysis of relay circuits. 
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ELECTRONIC PROGRAM-CONTROLLED COMPUTERS 


V.A. Los’ 
(Moscow) 


Digital methods have to be used in solving certain automatic control problems (e.g., those considered in [1}). 


In some cases it is best not to use digital computers however, because the computer is complex, and digit- 
alizing devices are required. Also, it is sometimes necessary to carry out many operations in order to solve a 
relatively simple problem. Combinations of digital and analog methods are therefore advantageous [2]. 


General-purpose and specialized analog computers have many units in common. For example, the MN-8 
has twelve multiplying units, about 40 nonlinear units, many operations amplifiers, etc. The desired type of 
machine should include a high-grade specimen of each of these types of unit. 


Switching (program) units enable one to use these units as required to solve equations. 


These computers differ from analog computers, although they have much in common with them, because 
they include digital devices. 


It is convenient to call computers of this type electronic program-controlled computers (EPC for short). 


They differ from analog machines in that they have all the features of digital machines, have many less 
operations amplifiers, and can be programmed to solve several distinct problems automatically, i.e., they can 
control several distinct processes at once. 


They can solve problems in logic, can select the optimum solution from many alternatives, and can solve 
problems that cannot be handled by analogs (e.g., can find the solution to a convolution integral). 


The design principles of such computers are presented here. 


1, EPC Working Principles 
Figure 1a illustrates the working principles of this type of machine. 


The computer is connected by the switches k,, ky, ..., kn, k', ky, . . . , ky to Solve the problem during 
certain time intervals t,. The program unit supplies to the computer,in the appropriate sequence voltage, pulses 
whose amplitudes are proportional to the inputs Uj,, Uj,, Ui, . - - , Uin- The pulses produced by the computer 
are suitably manipulated and are distributed,via the various channels, to the output units BY,, BY, . . . , BYn- 

The voltages reaching these units are stored, and are used in subsequent calculations or to control some 
process. 

Figure 1b shows (at 0) the pulses supplied by the time-interval generator TIG, which pulses control the 
program unit. Diagrams 1, 2, 3, ... , n show the outputs from the output units. The actual law followed by the 
function is shown by the broken line. Diagram £ shows the output from the computer unit. 


The principle is the pulse-height principle often used in time-division telemetry. The system differs from 
telemetry systems in that: 


a) the computer unit replaces the communication channel; 
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Fig. 1. 


b) the transducers, receivers and communication channel (computer unit) are combined into one; 


c) the transfer coefficient of the system is determined by the transfer coefficient o: the computer unit, 
whereas this coefficient must be constant if a communication channel is used, and so on [3, 4]. 


2. EPC Units 


Analog computer units are always present; their numbers are determined by the operations to be performed. 


These units include functional converters, differentiators, integrators, and so on. They are described in de- 
tail in a number of places, and so we shall not dwell upon them, 


We consider only certain special requirements on units used in EPC computers, and some particular aspects 
of the program-control units, 


a) The computing devices. Diagram © of Fig. 1 shows that these devices work as sampled-data units. 
However, any electronic system has a finite passband. Therefore, the frequency range of the controlled object 
is restricted by the fact that these devices are used repeatedly for different purposes. This frequency range is 
taken at present, so the passband and the pulse parameters must be chosen to correspond. 


The passband of such a device is controlled by the bandwidth of the wide-band operations amplifiers (OA). 
These amplifiers must be drift-free, to ensure stable working. The amplifiers designed at the Institute of Auto- 
mation and Remote Control [5] are adequate in this respect. The passband in the direct negative feedback mode 
is 10-15 kc for a signal of maximum amplitude 100 v. 
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| : b) The switches (gates), These are devices whose trans- 
fer coefficients k are near to one for control voltages of one 


Gy out sign (e.g., +Uy » —Uy,, Fig. 2), and zero for voltages of the 
opposite sign (—Uy,, +Uy,). By transfer coefficient we mean 
Ts, the ratio of the output voltage to the input voltage. 
Many types of electronic gates have been described. In 
Fig. 2. EPC we can use only gates whose linearity is good (whose 


transfer coefficients are constant) and independent of voltage 
and circuit parameters, and whose passbands are wide. The 
gates must not load the supply sources, and the control voltages 
demanded must not be large. 


A diode bridge with six diodes [6] satisfies this require- 
ment. 


The gate whose circuit is shown in Fig. 3 has been tested. 


The diodes are of silicon or hard-tube types. The 
+450-v lines must be stabilized. 


The tests gave these results: transfer coefficient 0.9265; 
deviation from the relation Uj.a), = 0.9265 Uj not more than 


Fig. 3. 0.05% over the entire range of inputs; passband up to 20 ke. 


Fig. 4. 


c) Program unit. This must supply the control voltages to the gates in the appropriate sequence. The out- 
put voltages must be of appropriate magnitudes, and of both signs. The stages must have power sufficient to con- 
trol several gates. 


A simple program unit is shown in Fig. 4; it uses the above diode gate. 


This unit contains a time-interval generator TIG, trigger circuits TC, a diode network DC, diode limiters 
DO, control triggers TY and cathode followers KII. The triggers are such as to provide plate voltage changes of 
#140 v, 


d) The output units. These units BY receive pulses of large mark-to-space ratio; 


Q = Tp/ty = 2NT/T = 2N, 


where T is the width of a control pulse, Tp is the time needed to solve the problem (1/Tp = fp = frequency with 
which the problem is solved) and N is the number of times the computer unit is used when the problem is solved 
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Fig. 5. 


These units store the pulses and deliver them up after a certain time for use in solving the problem. They 
also provide stepwise approximations to functions. Figure 5b and c show the approximated function f(t) and the 
approximating function f, (t), respectively. 


Each unit is an operations amplifier working as a delay line (when ky and kg of Fig. 5a are closed) and as a 
store (when ky and ke are open). 


3. Block Diagram of the EPC 
Consider how the system of equations given in [7] 


F, = asinz—bcosz = 0, 
Fy =acosz + bsinz —y =0, 


is solved; here a and b are inputs (time-functions), and x and y are the controlled variables. 
System (1) may, for convenience in use with computers, be put in the form 


F,=asinz—bcosz = 0, 


Fy=y—V =0. 


We do not need to follow the solution right through to the numerical results in order to illustrate the prin- 
ciples. No special symbols are needed to denote the scales in which the machine operates. 


Figure 6 shows the block diagram of an electronic device to 


ar M, sun solve system (1"). 
asing The symbols are: M, and M,, multipliers; sin and cos, cor- 
z ‘yas responding functional converters; 6, an amplifier of high gain, or a 
— Paty Js d difference amplifier; and 5, a summation point. 
b 
* = The device must be checked for steady-state stability when 
er}, the circuit is compiled [7]. 


Figure 7 shows a circuit to solve (1") by digital methods. 
The first equation of (1) is solved in two stages. 


Stage 1. The gates denoted by 1 are opened by the control voltages. The product,a sinx,,is fed to BY, 
where it is stored and inverted. 


Stage 2. The gates denoted by 2 are opened. The product, bcosx,, is fed to the difference amplifier 6. 
The stored product, asinx,, is also fed to this amplifier. 


The output x, from 6 varies in such a way that the input, asinxy + bcosx, is almost zero. If asinx, + 
+ bcosx, = 0, then x, is the first approximation to the solution of (1'). The root to the first equation in (1") is 
obtained by iteration. 
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Figure 8a gives a graphical interpretation of the way in which the first equation of (1') is solved. 


Let the device be in a state such that the output from 6 is x4, and let bcos x, be produced in the first stage, 
In the second stage, asinx, is compared with the voltage produced previously. The output from the amplifier is 
then a voltage x, such that asinx, — bcosx, is zero. This completes the first approximation. The process is then 
repeated, Several repetitions give us a voltage x,, that is the solution to the first equation of (1'). The second 
equation of (1°) is solved as a result of the solution of the first. 


No. of units in system 

No. of systems no 

of type (1) | program 
control 


OA 

Functional converters 
Gates 

Total diodes (roughly) 


OA 

Functional converters 
Gates 

Total diodes (roughly) 


OA 

Functional converters 
Gates 

Total diodes (roughly) 


a 
? 
2 
| 
| Fig. 7. 
asins -beosr asing -beosr 
aslar 
Fig. 8. 
| 
| 1 12 14 
6 9 
147 165 
| 2 15 28 
if 6 18 
20 - 
225 330 
| 3 18 56 
1 6 36 
36 
320 660 
| | aie 


The continuous-action system considered previously can be made 


p, U*% such that the machine always has available the solution, provided there 

Co is a solution that is unique in the mathematical sense. The program-~ 
“a max(a,b) controlled system may not always give the solution, although, formally 
D; speaking, one exists. Figure 8a shows an instance in which the solution 
> As 7} (a,b) would diverge. The solution diverges if the first stage produces a sinx, 

144 , and the second “ b cos xq. 

b : Some definite rules must be followed in designing such program~- 
+e controlled systems, For example, we must start to solve the first equa- 


tion in (1") by forming bcosx, and not asinx. 


Thus, the solution to (1") will converge only if the factors a and b ! 
are supplied in some suitable sequence. The lesser (in absolute value) 
must be the first to be supplied. 


A diode circuit (Fig. 9) may be used to separate the minimum and maximum factors. Of the pair Dy and 
D,, the one that receives the larger voltage will conduct; and of Dg and D,, the one that receives the lesser yolt- 
age [8]. 


It is possible to transform the equations in such a way that the number of iterations is minimized. 


Consider a graphical interpretation of one suitable method, as applied to the first equation of (1"). We 4 
transform the equation to 


x= 2z—asinz + bcosz. 


To solve this on the EPC we put it in the form 


+ asin 24, = bcos 
where i is the number of the iteration. 
Figure 10 shows that the rule given above must be obeyed in solving the transformed equation. 


In one instance, where a system of the type of (1') had to be solved, the program method required far less 
computer units than did a normal analog method. 


reasiar -(2+beosz) reasing -(2 +beosz) 
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Fig. 10. 


The table given below illustrates the gain to be obtained from program control when one, two or three sys- 
tems of the type of (1") are to be solved, Here we have assumed that each multiplier has three operations amp- 


i 
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lifiers and two squaring units, each of which contains 15 diodes, that the sin and cos units each contain two amp- 
lifiers andan arc~-sin converter (30 diodes), that the sin-to-cos converter contains two amplifiers and an arc-sin 
converter, that the gates each contain six diodes, and that the device for deriving moduli consists of three amp- 
lifiers and three squaring units. 
4. Some Technical Features of the EPC 

The following aspects require attention in designing such systems: 

a) the required repetition frequency (time) must be decided on; 

b) the signal-selection system at the input requires consideration; 


c) the necessary passband in the computer units must be calculated or, if this passband is fixed, the fre- 
quency range of the controlled object must be established. 


If we assume that the problem is solved in a single iteration, the period needed for a single solution is de- 
termined by the following factors: | 


a) the method of approximating the functions; 

b) the maximum permissible error in the approximation (AFg), and 
c) the way in which the functions and their derivatives vary. 
The error in the approximation is given by 


AF (t) = Fa(t) — f(t) = AFa, 


where F (t) is the function to be approximated and Fg (t) is the output produced by the output unit. 


Piecewise-linear and step-function approximations can be used. The first is the more accurate, but the ‘ 
function and its first derivative have both to be stored. The second is less accurate, but only the function has to 
be stored. y 
. If the second method is used, the time needed to solve the problem is 


Tp = 


The method of selecting the signals has to be considered, because, if the parameters of the program unit 
and of the operations amplifier are chosen wrongly, the transient from the previous pulse may not have died away 
when the next pulse arrives. The result will then be cross-talk between channels that should be separate. 1 


This selection problem can be solved as follows. We may assume that the mark-to-space ratio of the pulses 
is constant. 


The program control dealt with above gave an interval between pulses equal to the pulse length, t; (ty is 
the unused interval), i.e., the mark-to-space ratio is one. We may also assume (Fig. 1) that the pulse height & 
at the input to the computing unit will vary in a random fashion, Fj (t) = &jfj(t — 21T), where i= 1, 2,...,N. 
The signal then has a spectral width G; (jw) that is determined by the width of the spectrum from a single pulse [9}: 


x when 
= {6 when ¢>t, 


ot wt 
Gin (Jo) — hin (t) = at, 


The amplitude-phase response of the computer is Y (jw), so the output signal is a time-function 


fout (t)= \ Gin (0) ¥ (jo) 


—2 


(9): 


For simplicity we assume that the phase characteristic is zero, and that the frequency-amplitude charac- 
teristic may be represented by a step function; the expression can then be calculated from the series expansion. 


Knowing the output signal, we can make t, such as to reduce the cross-talk between channels to a suitable 

level, 
SUMMARY 

There are at present two ways in which analog computers are developing. 

The first is simply along the lines of electromechanical analogs. 


The second is new, and involves structural systems in which the computing units are used repeatedly to 
solve many equations with similar structures. Machines of this type are called electronic program-controlled 
analogs. 


This new type of machine economizes on computing units, and has the properties of a digital system; it 
can solve a far wider range of problems. 
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PHOTORESISTORS AS CIRCUIT ELEMENTS 


S.V. Svechnikov 
(Kiev) 


The main characteristics of commercial photoresistors are considered in relation 
to circuit uses in terms of the general laws of the internal photoelectric effect. It is 
shown that the best basic parameters to use are Ay = eyokBoTo/d’, and the coefficients 
y and a, which,respectively, define the nonlinearity of the voltage~current and inten- 
sity -current characteristics and the permissible power dissipation, and the mean carrier 
lifetime T». Methods of calculating the load and frequency responses, and temperature 
compensation circuits, are given. 


The equations enable one to calculate photoresistor circuits with an accuracy 
that is adequate for practical purposes. 


Developments in semiconductor technology in the last 5-10 years have made photoresistors normal elements 
in many signalling, switching, counting and automatic control circuits [1-16]. 


Photoresistors have much higher values of integral sensitivity and of permissible photocurrent than vacuum 
photocells do; they are also smaller. Their designs are simpler, and do not involve evacuated tubes. They work 
stably, and show no appreciable aging or drift. The average service lives of the usual commercial photoresistors 
are measured in tens of thousands of hours. 


The main disadvantages are the relatively slow response, the relation of response time to light intensity, 
the sensitivity of the dark current and the photocurrent to temperature, and the absence of saturation regions on 
the voltage-current curves, which latter makes the sensitivity dependent on the voltage applied. 


Table 1 and Fig. 1 give the basic parameters of commercial photoresistors. Data are also given in Table 1 
for photodiodes and small photomultipliers, for comparison. 


These data show that existing types are sensitive to the near ultraviolet, to visible light, and to near infra~ 
red radiation. Some types (FSK-M, FSS-M) are highly selective, others have broad sensitivity regions (FS-K, 
FS-D, FS~A, FS-B). Some [e.g., FSK-M, FSS-M, FS-K (copper doped)] are very sensitive toa-, B-, y~- and x- 


3 rays [1-4]. 


Table 1 also shows that photoresistors based on cadmium compounds have the largest integral sensitivity to 
visible light (e.g., FSK-M, FSS-M, FS-K, FS-D). These types can be made to have areas ranging from less than 
a mm to several cm?, and to carry currents ranging from microamps to amps. 


Photoresistors have several features that influence calculations on circuits in which they are used. These 
features have their source in the laws of the internal photoelectric effect. 


These laws are determined by the structure of the semiconductor, by any impurities present, by the nature 
of the optical transitions, and by the type of contact used. In general, the photocurrent Ip is related to the light 
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Fig. 1. Spectral sensitivity curves for: 1) FS-K1; 2) FS-K2; 3) FSK-M; 4) 
FS-B2; 5) FS-A1; 6) FSS-M; 17) FS-D. 


aA 
4 
Va 
3 4 5 6 7 InL, lux / 2 d 4 S5inU,v 
Fig. 2. Typical intensity -current curves Fig. 3. Typical voltage-current 
for photoresistors: 1) FS-K1, U = 50 v; 2) curves for photoresistors: 1) FS-K1, 
FS-K, U = 50 v; 3) FS-B2, U=40v; 4) L = 200 lux; 2) FS-K2, L = 200 lux; 
FSK-M, U = 20 y; 5) FS-Al, U= 15 y; 6) 3) FSK-M, L = 200 lux; 4) FS-B2, 
FSS-M, U= 13 y; 17) FS-D, U = 50 v. L = 200 lux; 5) FS-Al, L = 100 lux; 
6) FSS-M, L = 10 lux; 17) FS-D, L = 
= 200 lux. 


where Cy, Cy, By a*, B*, y* are constants for any given type of photoresistor. 


The experimental curves of Fig. 2 fit (1) very well; these curves are typical ones for commercial photore- 
sistors, Figure 2 shows that we can, for practical purposes, approximate a given intensity range by a power law 


Ip = AL®, @) 


where A and a are constants. 


Now I e expU for nonohmic contacts [5], and so we can [2] represent the yoltage-current relation by a 
power law also: 


Ip = BUY, (3) 


| Bi 
et 


where U is the voltage applied, and B and y are constants. If the contacts are ohmic, y = 1. Figure 3 illustrates 
(3); it gives curves for the usual photoresistors. 


The steady-state photocurrent in a semiconductor is [6, 7] related to L and U by 


eBk pt 4 


where e is the electronic charge, k is the absorption coefficient, 6 is the quantum yield, r and y are, respectively, 
the mean lifetime and mobility of the carriers, and d is the distance between the contacts, 


It has been shown (3, 6, 8] that, in general, r = r(L), 8 = B(L), w = w(U). The mobility is related to the 
voltage [3] by 


= (5) 
where p19 depends on the semiconductor and y' is a coefficient. 


We may relate the lifetime and quantum yield to L in a first approximation by expanding the functions in 
series form, assuming that the constant term is zero. We then have 


tL", B= BoL*’, (6) 
where a", a", To, 89 are constants. 


We substitute (5) and (6) into (4) to get 


Ip = AgL“UY, (7) 


Ao = (8) 


Now Ag, a and y are to be considered as parameters of the photoresistor; they can be used to calculate all 
the quantities given in Table 2. 


TABLE 2 
Parameters of Some Photoresistors 
Type 
anaimeter = 
4 <| 8 | gs-xe 
‘ Power dissipation, mw 10 | 100 | 100 150 10 | 400 | 30 
Voltage nonlinearity 
) coefficient, y 1 1 i i i* i i** 
Intensity nonlineari 
— 0.4} 0.5 0.8 | 0,5—0.8"**; 0.5] 0.8) 1.2 
A ximate carrier 
ifetime To, sec | 10-* | 10+ 10-* | 10-* | 10-* 
Temperature of 
total current, C, 
) for the range 20-100°C —0.6 1.2 0.8 0.3 0.8); 1.5] 1.5 


«Aluminum or cadmium contacts. 
**Cadmium contacts. 
***The lower value of a corresponds to lower intensities. 


| 
| 
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We see from comparing (7) with (2) and (3) that the coefficients of proportionality in the characteristics 
are, respectively, 


B= AgL*. (9) 


It is found that Ag is not strictly constant for photoresistors of any one type. The values of Ag may differ 
by 50% or more. The causes are variations in electrode spacing and chemical composition, particularly the latter. 
Equation (7) shows that Ag" 1/d*, Hence, Ag changes greatly with small changes in d. 


Types FS-K, FS-A and FS-D have relatively large values of d (3-5 mm), and the fluctuations are only 5 to 
10%, Type FSK-M (monocrystals) has d ~ 1 mm, and the electrode spacings may vary by 20-50%, 


>... 
Uy 
Fig. 4. Equivalent cir- Fig. 5. Derivation of the op- 
cuit of a photoresistor. timum load resistance. 


It is clear that d must be kept strictly constant if photoresistors are to be treated as interchangeable circuit 
elements. 


Besides Ag, we also have « and y as parameters; these specify the sensitivity and resistance. 
It is usual to give the static sensitivity, which is 


K,= (10) 


Now (7) and (10) show that 
K = Ag@~u, (11) 


Hence, K is not really a parameter, because the nonlinearity causes it to have a different value at each 
point on the curve. The mean K may be used in calculations on the permissible voltage for a given intensity, in 
which the permissible power rating Pmax is used. This power is found from the bend on the voltage-current curve 
(point P in Fig, 3). If P > Pygx, the temperature rises rapidly, the photocurrent falls, and the unit is damaged. 


Now (7) and (11) show that 
Pp max = [pUpmax = AoL"Upinax = KpLU pmax, (12) 
so 
Up max (13) 


Consider a small increase in photocurrent in relation to I, = f(L,U); we have 
Alp = p- (upAL + AUp), (14) 


where y = R,K,; Kj and Rj are the differential sensitivity and differential resistance, respectively, and are given by 


=const 


AU 


“4 R, = (a7), (15) 


Now K; and K are related by 
K/K =a, (16) 
a6 which formula should be used to calculate Kj. 
Figure 4 shows the equivalent circuit derived from (14). If 
a | R is the input resistance of an amplifier, the output parameter of 
the circuit in Fig. 4 is AU. Now 
Qs 
40 60 more AU = — AUy = — (17) 
Fig. 6. Relative change of current with 
the atL= eating (14) (27) we Reve 
= 100 lux and U = Uj, is taken as unity. “, 
1) FS-K1; 2) FS-K2; 3) FSK-M; 4) AU, = AL, (18) 
FS-M1; 5) FSS-M. 


where op, = Ry/ Rj. 
On the other hand, the load should receive the maximum power in the steady state. Figure 5 shows that 


as) 


where Rs = U/Ip is the steady-state resistance of the cell. Now dP};/dRy; = 0 occurs when R= Rg, 80 xzz = Ry/Rj. 
Then (3) and (15) give us that 


Re/Ri = 7, (20) 


which formula may be used to calculate Rj. Substitrting (20) into (18), we have 


Y 
The above argument shows that a photoresistor may be treated as a circuit element if we use the three basic 
parameters Ag, a and y. 
A photoresistor is, as a rule, the circuit element that shows the largest temperature coefficient. 


Photoresistor circuits have to be compensated against temperature changes at low light intensities. Figure 6 
shows the relative change of current with temperature for some standard types, The curves show that the tem- 
perature coefficient can usually be taken as roughly constant in the range from 0 to 100°C. 


Figure 5 then gives us that 


Ry (t°) = Rs (22) 


| 
by 


R,, (t,) 


(23) 
is given. 


Wetake R,, to consist of a resistance Ry of temperature coefficient ap in series with, say, a thermistor of re- 
sistance Ry that has a temperature coefficient ay. Then 


Ry (t°) = Ry (A + apAt®) + Ryo (1 


and Ry and Ry, are given by 
Ro + Rro = Ru (to); (24) 
Ro (i tarAt®) + Ryo (1 a,At°) = Ry (t°). (25) 
Therefore, 
R,(t,) (1 + a At*)—R, (0°) 
Ry = + (4p + (26) 
R, (ts) (1 + apAt®)— R, (2°) 
+ (4, + @p) (27) 
Now 8; = f(t*), so the relative error in the compensation is 
(t°)— 
(28) 


B 

Tests show that By is usually 1-2% over the range 20-100°C, The transient response, or the frequency re- 
sponse, is controlled by the mean carrier lifetime, rp. 

If the intensity-current curve is nonlinear, r will be a function of intensity. This follows directly from (2) 
and (4), Now (4), (6) and (7) give us that 

a’=a-~l1. 

Thus, if a < 1, we have a’ < 0 and a rise with intensity in accordance with r ~ T,/L™". Ifa > 1, we 
have a" > 0, and r falls as L rises. 

The formal kinetic equation for I,(t) can be put as 


4th (29) 


where 6* = ByekU/¢ is the quantum yield in current terms. 


The solution to (29) also gives us the frequency response. In general, r = r(L, Ip, t) and 8* = B(L, Ip, t). 
The kinetics will clearly then be nonlinear. The linear case occurs when r and 8* are constant. This, in fact, 
occurs only at high intensities, or if bias illumination is used. Any modulated beam has a certain mean illumina- 
tion level Ly. Sinusoidal [6] and square-wave [7] modulations have been considered; the formulas for Ip(t) are, 


respectively, 


Ip (30) 


| where 


and 
Ip = (77), 


which also give the frequency response. Now ry is the photoresistor parameter in (30) and (31). Table 2 gives 
the mean values of rT» for the standard types, and Fig. 7 gives the frequency-response curves recorded with square- 
wave modulation. The current Ino = 8* Tel corresponds to f=. 


| 10 100 1000 fcps 
Fig. 7. Frequency~-response curves for photoresistors: 
1 and 4) FS-K1; 2) FS-K2; 3) FSK-M; 5) FS~A1; 
6) FSS-M. 
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There is always a certain constant component in the illumination in such tests; this activates the specimen, 
and so we can assume that r = Tg = const. This activation effect will appear only at high values of L if single 
flashes are used. At low intensities we have 8* = 8*(t), as has been shown [4, 8,9]. The initial section of the 
current curve is then hyperbolic. In fact, (29) shows that, whent « 1, 


Tp (t) = B* (t) Lt. (32) 
We expand 8*(t) in series form, and get a general expression for the initial section of the intensity response: 
Tp (t) = BoLt”. (33) 
The special case 6* = const corresponds ton = 1. If n> 1, the photoresistor shows a delay in its response; 
this delay will be the larger the smaller L. 
This fact should be borne in mind in calculating photoresistor circuits, including correction circuits. 
) 
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’ THE STABILITY OF AN EXTRAPOLATION DEVICE 


Ludwig Prousa 


4 (Pardubice, Czechoslovakia) 


Much attention is at present being given to one type of sampled-data extrapolator, namely to the type that 
contains a sampling unit, summators, integrators and a feedback loop [1-10]. These devices are useful in very 
many different ways. 


This note deals with a special case (the one in which there are two integrators); it is, all the same, a very 
important case. It is the one dealt with in mest detail in the literature. 


1, Stability of an Extrapolator with a Gate 


Consider the device shown in Fig. 1 (the symbols are the same as Tsypkin's, so far as is possible). The gate 
has y « 1. The transfer function for the continuous-action part of the open-loop system is 
= 1 
a) 


We introduce the relative time, and put T,p = q, yaTy =a, ybT; =~ B, ts/T; = Tr; we get 


2) 
The original, corresponding to Y(q), is 


y(t) =0 when 


(3) 
y(t)=a+ when t>+, 


Thus we put e9 = z, and get the pulse-transfer function for the open-loop system as 


Y*(:)= (4) 
the characteristic equation is 
Br =0, (5) 


The stability conditions for the closed-loop system are 
B>0, 
a > (6) 
i 


»s|s 


Fig. 1. 


Thecharacteristic equation and the stability conditions are derived in another way in [5].* 


The special case r = 0 has been dealt with [9], where the stability region is split up into subregions that 
correspond to different forms for the roots of the characteristic equation. 


Theorem 1. Let the delay time r vary over the range 0 =< r < 1. Then the triangular stability region 
(that found for r = 0) is deformed (Fig. 2) in such a way that all the ordinates 8,on lines that split the triangle 
into subregions,remain unchanged for any r # 0; the abscissas a are 
displaced linearly by Br. 


The proof of this theorem follows from (5). Clearly, the point 
(a,8) for r = 0 and the point (a + 8T,8) for r = T correspond to 
identical roots. 


Theorem 2. The following inequalities give the bounds to a 
subregion, to which bounds correspond two real roots of the character- 
istic equation (one positive, one negative): 


B>0, 
a>i + Br, (7) 
). 


The proof of this theorem follows from (5). The necessary and 


sufficient condition for two real roots of opposite signs to exist is that the following inequalities are all complied 
with at once: 


[a + — > 4B, 
To = 
D4 of <0. 


It is readily seen that the second and third conditions in (8) are complied with at once only if 


«>1 + Br. (9) 


But the first condition in (8) is then also complied with. We add the first and third conditions in (6), and 
so get (7). 

Figure 2 shows three special cases, namely those for r = 0, r = Y2, r = 1 (see also [5]). The full lines cor- 
respond to the conditions of (6), and the broken lines to the second condition in (7). 


2. Stability of an Extrapolator with a Clamp 


In this case y = 1. We consider only the special case r = 0. Then 


*There are misprints in [5]; the term 8z in the characteristic equation has a minus sign, and the third stability 
condition should have < instead of >. 
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ie b 
We put aT, = a, bT; = 8, and get 
(= + 4). (11) 


Thus we put e9 = z, and get the pulse-transfer function for the open-loop system as 


i i 
Y*(z)= [e+ (11) 


Theorem 3. The circuit with a clamp with r = 0 is equivalent to the one with the gate when r = 5. 


+ 1-6 4? 


Fig. 3. 


The proof of this follows from (12) and (4) taken with r = 44. 
The clamp circuit (Fig. 3) has been considered in detail [2,10]. In [10] use is made of the parameters 


and the bifurcation of the stability region is given for these. Now a and 86 are related to j and k: 


(14) 
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A PEAK-HOLDING REGULATOR THAT FOLLOWS THE PEAK 


N.V. Grishko 
(Moscow) 


A peak-holding regulator that follows the peak is described. It consists of two 
main parts, one of which measures the distance from the peak, and the other of which 
is a sampled-data follow-up device. Test results are given; these tests had as object 


to find parameters optimal in the sense that the root-mean-square follow-up error is 
minimal. 


The extremum (maximum or minimum) has to be sought in cases in which the peak changes its position 
continuously. The response of the system is improved if a servo is used to follow the extremum. 


Fig. 1. The experimental system for studying a peak-holding regulator that 
follows the peak. 


The author's system, which is described here, used the signal that indicates the deviation from the peak as 
the error signal for the sampled-data servo, A regulator of this type is under development for industrial purposes; 
a pneumatic system is to be used. An electronic analog was first tested.* The controlled object and the device 
for adjusting the regulator to the peak were included in the analog. The tests showed that the system had some 
advantages over other peak-holding devices that have been described. 


*From proceedings of the experimental seminar on pneumo- and hydroautomation control, IAT, AN SSSR, May 
22, 1958. 
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1. The electronic regulator, This (Fig. 1, BA) has three main 
parts, namely the mismatch measurer (A), the sampled-data servo 
@) and the control unit (not shown). 


Unit A consists of an oscillation source (A1), the measurer it- 
self (A2), and a floating-zero read-off system (A3). 


Unit Al of Fig. 1 supplies a train of square pulses to the input 
of the controlled object. The pulse length and train length can be 
varied, The output signal produced from the object by these pulses 
gives the magnitude and sign of the derivative of the steady ~state 
characteristic for the point at which the object actually is. This 
output is fed to two integrators (4 and 10 in Fig. 1). One integrator 
deals with the output when positive inputs are supplied to the object, 
the other with the output corresponding to negative inputs, The out- 
puts from the integrators are differenced, for which purpose an in- 
verter and summator are used (amplifiers 11 and 6 in Fig. 1). The 


Fig. 2. Diagram of the working se- resulting output signal (Fig. 1,A and Fig. 3, curve 4) specifies the 
quence of the contacts in the control distance from the peak. This method, as is usual, reduces the effects 
unit, Ty, is the time taken to meas- of noise. The fixed component in the output from the object is 

ure the distance from the peak, Tp is balanced out up to the start of the measurement by connecting in 
the pulse length in the servo, Tf is the floating system by closing contact ky. At the same time, con- 
the time for which the floating-zero tacts ky and ky are closed to discharge the integrators 4 and 10. 


system operates, and during which the 
integrators in A2 discharge, and T, is 
the time for one working cycle. 


The resulting output signal is fed to the servo by closing kg. 
The servo drives the controller in two distinct ways; it causes it to 
advance by steps, and to drive at a constant rate between steps. The 
sizes of the steps in position or rate depend on the error signal fed to 
the servo. 


The system works in the above way when contacts ky, ... , kp are closed in accordance with Fig. 2, which 
' describes the operations of the control unit in the regulator, which here consists of the control unit of the elec- 
tronic analog and of a set of contactors. 


Fig. 3. 1) Output from object; 2) position of controller corresponding to the extreme 
value of the output; 3) position of controlier during follow-up; 4) distance from peak; 
5 and 6) outputs from integrators 4 and 10. 
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Fig. 4. 1) Root-mean-square follow-up error 
when the peak is moving and noise is present; 
2) the same, but without noise; 3) mean- 
square value of the pulses fed to the sampled- 
data servo when the peak is moving rapidly 
and noise is present; 4) the same, but without 
noise; 5 and 6) correspond to 3 and 4, but with 
the peak moving slowly. 


2. Adjusting the regulator. The parameters of the — 
pulse train can be adjusted; so can the position and rate 
steps produced by the servo in the controller. These ad- 
justments are made manually to give minimum follow-up 
error under given conditions. The optimum parameters 
were selected by measuring the root-mean-square (integral) 


error. Unit B’ in Fig. 1 was used for this purpose. The ac- 


tual position of the peak is, of course, unknown under real 
conditions; therefore, unit E (Fig. 1) measured the mean 
value of the square of the pulse input to the servo, which is 
a quantity that can be found in practice. There would, in 
principle, be no difficulty in making the adjustments auto- 
matically. The peak-holding regulator could be adjusted 
continuously by an analogous peak-holding regulator, which 
would seek out those values of the parameters for first reg- 
ulator that result in the minimum mean-square input to 

the servo. 


3, Results. The analog to the object (Fig. 1,C) was 
a squaring unit with lags before and after it (C1), a ball 
noise generator (C2) and a source (C3) that produced saw~ 
tooth oscillations in the position of the peak without affect- 
ing the peak height. The parameters of these various units 
could. be altered. Figure 3a shows the results for zero noise. 
Curve 3 shows the motion of the controller; the train of 
pulses (one), the step, and the period of operation of the 
floating system are clearly seen. The left section shows 
what happens when there is no continuous motion between 
the steps. If we compare the left and right halves of curve 1 
in Fig. 3a, we see that the continuous motion of the con- 
troller improves the response; the follow-up worsens around 
the extrema on curve 2 of Fig. 3, which shows that the 
servo system stands in need of improvement. Figure 3b 


and c show the follow-up when noise is present; the filtering action of the error measurer may be judged from 
these curves. The relation of root-mean-square error to the parameters was calculated from the results. Figure 4 
shows an example (the parameter that was varied here was the number of pulses in a train); there is a marked op- 
timum point which moves to the right as the noise levei increases or the rate of movement falls. 


I wish to thank T.N. Kolerova and M.M. Sharafetdinoy for assistance in the experiments. 


Received October 24, 1958 


METHODS OF STUDYING THE STABILITY OF AUTOMATIC 
CONTROL SYSTEMS WITH DISTRIBUTED PARAMETERS, 
WITH ALLOWANCE FOR LOSSES 


Ya.B. Kadymov 
(Baku) 


The stability criterion for systems with distributed parameters [1] is extended to 


systems with losses. 


Consider automatic control systems with distributed parameters whose characteristic equations can be put 


in the form 


sve 


av 
* =0, 
where #;(p), ¥:(p) are polynomials of degrees n and m, with m <n, 
Y=V + pb, 


or, in the general case, 


y=V ph+a, 


where 1, c, b and a are the parameters of the part with distributed constants. 
Consider the function 


where zeros equal the roots of (1), and whose poles equal the roots of 


¥1(p) = 0. 


We use the notation 


(1) 


(2) 


(3) 


(4) 


(5) 


(8) 


di(p)e* — va (p)e _ Sal?) 
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It is readily shown that ¢,(p) is an analytic function. Then the critical values of Ig and wei can (2) be 
found from 


(ie) = 1 g (8) 
or 


Let us transform the expression for yl: 


Here r = 1//c is the wave propagation time. 
Now A (p) can be put in the form 


(p) = tp V 145. (11) 


We expand /(I + b/cp) as a series [3] and get 


ib i 4 ix3 1 
The series converges if |1/p| = 1. This condition is obeyed with p= jw if |w| > 1. 
Thus, if |w| > 1, we have 


But A (p) can be put in another form 


A(p) ex V 1+ 
We expand 7 (1 + cp/b) and get 


This series converges if |p| <1, i.e., with p = jw the expansion applies if |w| < 1. Now #1 = 2/2(1+)), 
so if |w| < 1, we have 


We need take only the first two terms in (13) and (15). If |w| > 1, we use the following expression for X (jw): 
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Then 

|| 
) 
) 

) 


(jo) = «(jo + (16) 


Then (11) and (16) give 
where 
ib 
m= 
Hence W (jw) can be put as 
W (je) = | W (ja) |e) Ke) (18) 


where | W (jw)| is the modulus of W (j,w), and arc tg f(w) is the argument of W (j,w). 
We substitute (18) into (17) and get 


|W (joo) | jor, (19) 
Taking logarithms of both parts of (19), we have 
In| W (ja) | + j arc tg f (@) = — mt — jor. (20) 


Separating the real and imaginary parts, we have 


In| W (jo)| = —mr, (21) 
or 
In| (23) 
——arctg f(a). (24) 
If |w| < 1, we will have 
(ja) = (1 +) (1+ few). (25) 


We perform the same expansions as for | w| > 1, and get for |w| < 1 that 


4m 
= (27) 


We give w values from —oo to a, and construct the r = f,(w) curve from (23) and (26), and the r = {(w) 
curve from (24) and (27). The points where these curves meet define w,,; and T,,;. 


We define the stability region from the limiting characteristic (for r = 0) in the way given in [1, 2]. 
We now demonstrate how to find w,,; and Tq Without resort to expanding A(p) = yl.in series form. 


») 


1) 
2) 


3) 


4) 


5) 


1) 


We put @* — jo = in the form 


We replace (11) by 


—io.8 arc tg —- 
|W (jo) | ef are te f ot cos 


Taking logarithms of both parts of (29), we have 


In | W (je) | + arctg f = — jt (2) 


wt + ( {in | W (je) | + jare tg f (@)}. 


We give w values from —oo to oo, and construct the curve 


b 
(ja) | W (je) | + fare tg f (w)). 


The points where (jw) meets the imaginary axis define wWorj and Teri. 
The regions of stability and instability are then derived as above. 


(29) 


(31) 


(32) 


The same method can be used to study the stability of systems that have (3) as their expressions for y. 
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